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SECTION  1 


INTRODUCTION  j 

I 

j 

j 

This  report  is  the  second  interim  technical  report  under  the  { 

Adaptive  Antenna  Control  program, Contract  No. DAAB07-76-8-85.  The 
Adaptive  Antenna  Control  (AAC)  prograun  is  under  the  contract  manage- 
ment of  the  Project  Manager  for  DCS  (Army)  Sojnmunications  Systems, 

Ft. Monmouth, New  Jersey.  E^OM, Ft. Monmouth, N.J.  is  responsible  for 
technical  direction  of  the  program.  SIGNATRON, Inc.of  Lexington, 
Massachusetts  with  RF  Systems,  Inc.  of  Cohassett  as  a major  sub- 
contractor have  the  responsibility  for  executing  the  program 
tasks . I 

The  AAC  program  encompasses  the  investigation  of  adaiptive 
antenna  control  techniques  to  enhance  communication  on  strategic  j 

trans-horizon  radio  paths.  For  paths  dominated  by  diffraction  | 

effects  with  a small  or  non  existent  scatter  component,  the 
emphasis  is  on  antenna  steering  techniques  to  mitigate  beam 
blockage  due  to  changes  in  the  refractive  index  gradient.  Anal- 
ysis of  the  diffraction  problem  has  been  included  in  the  first 
Interim  Technical  Report  Ll.l]  and  a design  of  a mechanically 
steered  feedhorn  with  electronic  adaptation  was  presented  in  the 
AAC  Design  Plan  [1.2].  The  implementation  and  field  test  of 
this  system  is  not  part  of  the  existing  contract. 

For  trans-horizon  paths  dominated  by  forward  scatter  recep- 
tion of  conmiunication  signals,  the  failure  to  successfully  il- 
luminate the  useful  scattering  volume  results  in  an  aperture  to 
medium  coupling  loss  and  neglects  an  important  diversity  capabil- 
ity associated  with  signal  angle  arrivals  from  different  portions 
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of  the  scattering  volxime.  This  program  includes  a comprehensive 
modeling  of  this  phenomenon,  development  of  a prediction  model, 
analysis  of  angle  diversity  capability,  design  and  development 
of  a combiner  system  for  augmentation  to  8th  order  diversity 
using  angle  of  arrival  signals,  design  and  development  of  a dual 
vertical  feedhorn  modification  for  the  C-band  RADC  test  system, 
and  an  extensive  angle  diversity  test  program  to  be  accon^lished 
on  the  RADC  Youngs town- Verona  TRC-132  test  link  in  the  fall  and 
winter  of  1977.  The  first  Interim  Technical  Report  [l.l]  de- 
veloped much  of  the  theoretical  basis  for  the  troposcatter  model 
and  the  prediction  technique.  The  Design  Plan  [1.2]  contains 
details  of  the  combiner  and  feedhorn  systems,  and  the  Test  Plan 
[1.3]  describes  the  factory  and  field  tests  in  detail. 

This  report  continues  with  the  theoretical  development  of 
the  troposcatter  model  and  prediction  technique.  Section  2 re- 
views other  approaches  to  troposcatter  modeling  and  develops  the 
basis  for  predicting  median  path  loss,  2(7  multipath  spread, 
aperture  to  medium  coupling  loss,  and  signal  correlation  distances 
at  the  receiver  antenna.  Angle  diversity  systems  are  considered 
in  Section  3 where  the  methods  of  Section  2 are  used  to  derive 
system  design  parauneters.  The  variation  of  hourly  median  sig- 
nals is  also  treated  in  Section  3 and  the  approach  for  determin- 
ing system  availability  is  outlined.  The  performamce  of  angle 
diversity  systems  with  digital  signaling  is  treated  in  Section  4. 
A method  of  computing  the  average  bit  error  probed>ility  for  a 
troposcatter  adaptive  equalizer  in  the  presence  of  both  in^licit 
and  angle  diversity  is  derived  for  the  first  time.  This  method 
is  used  in  Section  4.2  to  confute  system  performance  for  some 
example  links. 
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Some  of  the  major  accon^lishments  and  findings  developed  in 
I this  study  zure 

* aperture  to  medium  coupling  loss  is  best  considered  an 
integral  pzurt  of  the  path  loss  rather  than  a separable 
entity  to  be  calculated  by  itself, 

• a numerical  integration  technique  derived  from  tiirbulent 
scattering  theory  has  been  developed  to  confute  the  path 
loss  including  aperture  to  medium  coupling  loss  effects, 

• an  accurate  closed  form  approximation  of  the  above 
integration  technique  has  been  developed  for  the  wide- 
beam  path  loss, 

• a closed  form  relationship  for  aperture  to  medium  coupl- 
ing loss  has  been  developed*  TtiB  widd»eam  path  loss 
and  aperature  to  medixim  coupling  loss  equations  when 
used  together  show  good  agreement  with  the  numerical 
integration  path  loss  result, 

* vertical  angle  diversity  is  superior  to  horizontal  angle 
diversity  for  typical  narrow  beamwidth  strategic  tropo- 
scatter  systems, 

* the  optimum  squint  angle  of  a duplex  dual  vertical  angle 
diversity  system  is  approximately  one  beamwidth.  The 
optimum  boresight/horizon  euicle  is  about  1/4  beamwidth 
for  L-band  systems  and  1/2  beamwidth  for  C-bzmd  sys^ 
terns. 


there  is  a long  term  advantage  for  euigle  diversity  over 
space  diversity  due  to  the  decorrelation  of  the  scat- 
tering volumes  in  the  former.  This  advantage  is  reduced 
by  the  greater  variance  of  signal  strength  from  the 
elevated  beam.  Bi^irical  data  from  a Bell  laboratories 
study  [1.4]  Indicate  that  the  two  effects  tend  to  cancel 
each  other  out. 


a method  requiring  the  calculation  of  a determinant  has 
been  developed  for  the  computation  of  average  bit  error 


prbbability  of  a troposcatter  equalizer  in  an  angle 
diversity  system  with  multipath  dispersion, 

* on  a system  comparison  basis,  dual  space/dual  frequency 
(2S/2F)  is  generally  inferior  to  dual  space/dual  angle 
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(2S/2A) . The  superiority  of  2S/2A  results  chiefly  from 
the  3 dB  advantage  of  the  second  power  amplifier  exceed- 
ing half  the  sguint  loss  in  dB  of  the  elevated  beam*  The 
implication  is  that  the  number  of  frequencies  required 
for  the  strategic  troposcatter  plant  could  be  halved 
and  overall  system  performance  would  be  improved, 

• the  preliminary  results  on  aperture  to  medium  coupling 
loss  as  a function  of  frequency  suggest  that  system 
improvement  may  result  if  L-band  systems  are  converted 
to  S-band. 

The  major  portion  of  the  remaining  study  effort  will  con- 
centrate on  validating  the  median  path  loss  predictions  in 
order  to  provide  an  accurate  absolute  median  path  loss  measure 
and  to  verify  the  coupling  loss  dependence  with  frequency.  The 
results  of  the  field  tests  will  also  provide  an  important  em- 
pirical guide  for  path  loss  prediction  and  for  evaluating  the 
long  term  decorrelation  advantage  of  vertical  angle  diversity. 
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SECTION  2 


TROPOSCATTER  PERFORMANCE  PREDICTION 

2.1  Introduction 

In  this  section  we  shall  describe  in  detail  the  background, 
theoretical  as  well  as  historical,  for  the  development  of  the 
SIGNATRON  trc^oscatter  software  package.  This  discussion  will  include 
a theoretical  development  of  the  propagation  model, analytical  per- 
formance prediction,  and  some  numerical  results.  A comparison 
with  the  results  of  other  authors  will  also  be  presented. 

Theory  and  practical  techniques  for  troposcatter  communica- 
tion started  developing  around  1950.  [Booker  and  Gordon,  1950.] 
Two  main  theories  have  been  proposed,  the  layer  reflection  theory 
[Friis,  et.  al.,  1957]  and  the  turbulent  scattering  theory. 

The  turbulent  scattering  theory  is  based  on  the  theory  of  tur- 
bulence structure  developed  in  1941  by  Kolmogorov  and  Obukhov. 

The  most  complete  treatment  may  be  found  in  Tatarski  [1971],  The 
layer  reflection  theory  received  the  most  attention  in  the  beginning, 
partly  because  it  predicted  the  observed  linear  dependence  of  the 
scattering  crossection  with  wavelength.  This  theory  also  predicts  a 
dependence  on  the  scattering  angle  of  the  form  e“^.  In  later 
years , experiments  at  higher  frequencies,  and  hence  narrower 

antenna  beams,  have  confirmed  the  validity  of  the  turbulence 

-11/3 

theory,  with  a scattering  angle  dependence  of  the  form  e 

and  a wavelength  dependence  for  the  scattering  crossection  of 
-1/3 

the  form  \ . In  practice  both  phenomena  can  exist  simul- 

taneously with  layer  reflection  being  more  frequent  when  wide- 
beam  antennas  and/or  low  frequencies  are  used.  Experimental 

results  fEklund  and  Wickerts,  19681  indicate  that  the  wavelength 

-1  3 

dependence  cah  vary  between  \ to  x . The  scattering  angle 

*2  —6 

dependence  has  been  measured  to  lie  in  the  range  of  e to  e . 
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The  exponent  in  the  scattering  angle  is  called  the  refractive 
index  spectrum  slope  since  it  also  indicates  the  fall-off  of  the 
wave  number  spectrum  of  the  refractive  index,  provided  the  wave- 
length is  in  the  so-called  inertial  subrange  [TatarsTci,  197l,p52'1. 

A third  theory  that  has  been  proposed  is  the  normal  mode 
theory  [Carroll  and  Ring,  1955%  This  technique  consists  of 
considering  the  atmosphere  as  a waveguide  propagating  infinitely 
many  modes  of  the  electromagnetic  field.  Except  for  diffraction 
path  analysis  this  method  is  generally  not  practical  to  use. 

A number  of  theoretical  and  computational  models  have  devel- 
oped [Rice,  et  al,  1967;  Hartman  and  Wilkerson,  1959;  Booker  and 
Gordon,  1950;  Tatarski,  1971;  Yeh,  19603,  Further  references  and 
a discussion  of  these  models  can  be  found  in  Larsen  (1968)  or 
Ranter  (1972) . Most  of  these  models  have  been  constructed  to 
agree  with  empirical  measurements  of  a few  selected  parameters, 
parttularly  the  received  signal  level  and  its  daily  or  seasonal 
variations.  As  a result  the  majority  of  the  models  are  not 
suitable  for  prediction  of  a nuniber  of  other  parameters  of  interest, 
such  as 

* coupling  loss, 

* delay  spread, 

* receiver  correlation  distances* 

* fade  rate* 

* angle  diversity  correlation* 

* the  dependence  on  distance, 

* elevation  angles, 

* frequency  , and 

* standard  atmospheric  parameters. 
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The  SIGNhTRC^  troposcatter  computer  program  is  general  in  scope 
and  can  handle  most  of  the  parameter  of  interest  without  nudcing 
unreasonable  assumptions  about  the  antenna  gain  patterns.  This 
greater  generality  also  means  that  the  path  has  to  be  specified 
in  a way  more  consistent  with  the  underlying  physics  than  has 
typically  been  done  in  the  past.  The  ensuing  discussion  will 
therefore  center  on  establishing  a firm  basis  for  describing 
the  troposcatter  model  and  its  associated  performance.  This  in- 
cludes separate  treatments  of  the 

• properties  of  the  troposphere 

• path  geometry  and  antenna  patterns 

• walytical  evaluation  of  the  model 

• assessment  of  the  effects  of  the  approximations  made 

• cougar ison  with  analytical  and  empirical  results. 

The  majority  of  the  discussion  will  be  directed  to 

• small  to  medium  distance  path 

• microwave  frequencies  (300  MHz  to  20  GHz) 

• horizontally  stratified  atmosphere. 

Correction  terms  for  long  distance  paths  are  included.  Most  of 
the  results  apply  to  higher  frequencies  than  listed  above,  but 
multiple  scattering  and  large  atmospheric  absorption  will  limit 
the  direct  application  of  the  model  at  higher  frequencies. 

2.2  The  Troposphere  for  Bevond-the«HOrlzon  Propagation 

At  the  frequencies  of  interest  the  tropospheric  radio  pro- 
pagation is  well  described  by  the  theory  of  geometrical  optics. 
We  shall  thus  primarily  use  the  ray-propagation  technique  and 
will  ignore  those  almost  line-of-sight  paths  vdiere  diffraction 
is  the  main  mode  of  communication. 
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This  section  is  intended  partly  as  a reference  for  the  sub- 
sequent sections  and  partly  as  an  introduction  to  some  of  the 
theories  of  atmospheric  structure  employed  in  tropospheric  scatter 
performance  prediction.  The  topics  to  be  discussed  in  this  sec- 
tion are  the  standard  atmosphere  and  refractive  index  variations 
in  Section  2.2.1,  layer  reflection  in  Section  2.2.2,  refractive 
index  fluctuations  in  Section  2.2.3  and  turbulent  scatter  in 
Section  2.2.4. 

2.2.1  The  Standard  Atmosphere 

Radio  meteorological  measurements  normally  yield  pressure, 
temperature  and  water  vapor  content  of  the  atmosphere.  From 
such  measurements  the  refractive  index  n,  or  the  refractivity 
N *=  (n  - 1)  • 10^,  can  be  determined  from  an  empirical  formula. 

The  relation  most  commonly  used  is  rsmith  and  Weintraub,  1953], 

N - (n-1)  10®  = (P  + 4810  • ^),  (2.1) 

where 

p * total  pressure  in  millibars 
e * water  vapor  pressure  in  millibars 
T • temperature  in  Kelvin  degrees. 

This  foxmula  is  designed  to  be  valid  in  the  temperature  range 
-50°  C to  40°  C with  an  error  less  than  .5%  in  N up  to  30  GHz. 

An  inqproved  equation  has  been  suggested  by  Thayer  (1974)  but  for 
most  radio  application  Eq.  (2.1)  is  sufficiently  accurate. 


r 


1 


Measurements  of  the  refractivity  as  a function  of  height 
above  the  surface  of  the  earth  results  in  a typical  profile,  as 
shown  in  Fig.  2.1.  The  refractivity  at  the  surface  is  of  the 
order  of  300,  and  a nearly  linear  decrease  with  height  is  normally 
found  in  the  first  two  or  three  kilometers.  In  the  standard 
atmosphere  the  slope  is  40  N-units/km.  This  line  is  also  sho«m 
in  Fig.  2.1. 

The  nearly  linear  decrease  with  heights  at  low  altitudes  makes 
it  convenient  to  introduce  an  "effective  earth  radius"  to  account 
for  the  ray  bending.  We  have  shown  previously  [Monsen  et  al,  1976] 
that  modeling  the  refractive  index  in  the  form 

r Rq  Y 

n(R^  + h)  « ^(Rq)  L ^ ^ J » o < h < 2 km 

o 

leads  to  the  exact  relationship 


In  these  equations  « 6368  km  is  the  actual  earth  radius,  h 
is  height  above  the  surface  and  y a parameter  determining  the 
slope  of  the  refractivity  profile.  R^  defines  an  equivalent  earth 
radius  in  a transformed  coordinate  system  with  preserved  distance 
measure  on  the  surface  and  with  straight  line  propagation  in  the 
atmosphere.  For  the  standard  atmosphere  y is 


y [standard  atmosphere]  * 0.255, 


corresponding  to  the  earth  radius  factor  K given  by 


e 4 

K • ~ ~ “ for  the  standard  atmosphere, 
o 
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CCIR  even  defines  a standard  refractivity  profile  of  the  form 


N 289  e 


-0.136  h 


In  this  form  N(h)  is  close  to  the  observed  profile, but  for  most 
troposcatter  applications  it  is  adequate  to  use  the  linear 
approximation  valid  near  the  surface.  The  height  of  the  conmon 
volume  on  a troposcatter  path  of  length  d vihere  the  antennas  are 
aimed  at  the  horizon  is 


h ~ 
c 


2 Ion  if  d B 370  1cm  (230  mi) 
and  K “ 4/3. 


For  longer  paths  an  equivalent  value  of  K may  be  found  so  that  the 
scattering  angle  is  correct.  The  small  altitude  errors  introduced 
by  this  approach  can  be  ccxnpensated  for  in  the  final  analysis. 

During  conditions  of  super  refraction  it  can  sometimes  be  pos- 
sible to  have  a direct  path  between  transmitter  and  receiver  on  a 
path  where  trcsposcatter  would  normally  be  the  only  mode  of  coinnuni- 
cation.  This  can  occnir  when  the  gradient  of  refractivity  is 
steeper  than  -157  N-units/km.  In  generaL  variations  of  surface 
refractivity  and  gradient  are  functions  of  time-of-year  and  geo- 
graphical location.  Data  have  been  accunmiulated  over  many  years 
to  allow  a reasonable  prediction  of  these  parameters. 

2.2.2  Reflection  frcxn  Atmospheric  Lavers 

The  predcsminantly  horizontal  stratification  of  the  atmosphere 
often  contains  local  steep  gradients  in  the  transition  frcxn  one 


value  of  N to  another.  For  electromagnetic  waves  with  a wavelength 
large  compared  to  the  width  of  the  layer  boundary  such  high 
gradient  regions  look  like  an  abrupt  transition  and  can  cause  a 
substantial  reflection  of  energy.  At  higher  frequencies  the 
transition  appears  more  gradual  and  little  or  no  energy  is  reflected. 
At  intermediate  frequencies  layer  reflection  and  turbulent  scattering 
can  exist  simultaneously,  greatly  complicating  the  analysis  of 
troposcatter  links.  In  this  section  we  briefly  summarize  some 
results  relating  to  layer  reflection  or  "feuillet*'  scattering. 

Layer  reflection  and  scattering  were  analyzed  by  Friis,  Crawford 
and  Hogg  (1957).  The  effective  scattering  crossection  was  found 
for  several  types  of  layer  structures.  The  dependence  of  the 
crossection  a^,  on  wavelength  X and  scattering  angle  is  listed 
below  (from  Rice,  et.al.,  1967): 


Large  Layers: 

6-" 

(2.2) 

Intermediate  Layers: 

x^ 

e"* 

(2.3) 

Small  Layers: 

x" 

6-“ 

(2.4) 

This  may  be  compared  with  the  mathematical  formula  listed  below, 
based  on  empirical  results  from  radio  data  [Norton,  1960;  Rice 
et  al,  1967], 

-5 

a ■ X 0 • 5“  » (2.5) 

o 

2 

where  a is  the  variance  of  the  refractive  index  and  L is  the 
n o 

outer  scale  of  turbulence  (see  Section  2.2.3  for  further  details). 
This  formula  is  close  to  the  behavior  expected  from  intermediate 
size  layers.  More  recent  high  frequency  data  (>16Hz)  indicate  a same- 
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what  different  form  of  a $ closer  to  that  predicted  by  the 

s 

Kolmogorov- Ctoukhov  turbulence  theory  described  in  the  next  section 


r 

\ 


L 


2.2.3  Refractive  Index  Fluctuations 

2. 2. 3.1  The  Covariance  Function 

Superimposed  on  the  average  refractive  ind^  variations  dis- 
cussed in  Section  2.2.1  is  a small,  rapidly  fluctuating  term, 
attributed  to  turbulence.  Let  n(r)  be  the  refractive  index  at  a 
point  r,  and  with 

n(r)  “ n (r)  + n (r) 

— o — 1 

where  n^(r)  ■ E (n(r)}  and  n^^  (r)  is  the  small, 

fluctuating  con^onent, 

n(£)  can  be  considered  as  a random  field.  The  covariance  function 
<Pn 


cp^(r,  r')  * E {n  (r)  n (r')T 
n — 1 — 1 


If  cp^  is  a function  only  of  r - r'  then  the  random  field  of  n(r) 
is  said  to  be  homogeneous , and  the  covariance  function  is  written 
as  - r').  If,  in  addition,  is  only  a function  of  Ijc  - jc' | 

then  the  random  field  is  said  to  be  isotropic  (in  addition  to 
homogeneous ) . The  wavenuniber  spectrum  of  a homogeneous  field  is 
defined  as  the  Fourier  transform  of  the  covariance  function. 


’n<2£> 


’'2-  ’'3> 


(2lt) 


JJJ 


i X*r  .3 
e dr 


2. 2. 3. 2 The  Structure  Function 

^e  structure  function  D (r,  r')  is  defined  by 

n *“  ““ 
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D^(r,  r'J  » E ■[[nj^(r)  - n^(£’)]^]-  . 

Where  £' ) is  a function  of  jc  - r'only,  the  random  field 

n(r)  is  said  to  be  locally  homogeneous . If  £' ) is  a func- 

tion of  Ijr  - X*  I only,  then  the  field  is  locally  homogeneous  and 
isotropic.  The  structure  function  is  occasionally  more  conyenient 
to  use  than  the  coyariance  function  since  it  allows  a simple 
description  of  locally  homogeneous  fields  without  assuming  homo- 
geneity. Boweyer,  the  atmosphere  is  neither  homogeneous  nor 
locally  homogeneous  in  the  sense  aboye.  At  a giyen  point  in  space 
it  is  instead  possible  to  find  a small  surrounding  yolume  within 
which  the  condition  for  homogeneity  is  satisfied,  and  a slightly 
larger  volume  within  which  the  condition  for  local  homogeneity  is 
satisfied.  In  most  practical  cases  the  extra  generality  of  the 
structure  function  is  not  very  important.  In  a nearly  homogeneous 
yolume  the  two  functions  are  related  by 

D (r)  - 2<p  (0)  - 2cp  (r),  (2.6) 

n n n 

where  _r  n^  represents  the  difference  r - _r'. 

2. 2. 3. 3 Scales  and  Structure  Constant 

In  terms  of  the  aboye  defined  functions  we  can  define  two  important 

scales  of  the  turbulence.  The  inner  scale  l , of  the  turbulence 

____  o 

is  the  dimension  of  a small  region  near  £ ^ o where  the  structure 
function  is  well  approximated  by  the  first  term  in  its  Taylor 
series.  If  the  turbulence  is  isotropic  this  scale  is  the  same 
in  all  directions.  The  outer  scale,  L , is  the  correlation  dis- 

— — — Q 

tance  of  cp^(r).  These  scales  will  be  tied  to  specific  formulas 
later,  allowing  us  to  eyaluate  them  numerically.  In  isotropic 
turbulence  with  £ in  the  range 

-t  « r « L 
o o 
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we  can  usually  write 


D (r) 
n 


2 q 
C 

n 


(2.7) 


The  constant  C is  called  the  structure  constant  of  the  refractive 
n 

index,  and  plays  a very  important  role  in  the  study  of  turbulence. 
The  power,  q,  has  been  shown  by  the  theory  of  Kolmogorov  and 
Obukhov  [Tatarski,  1971]  to  be 

q - 2/3. 

2 .2.3.4  The  von  Karman  Spectrum  and  the  Kolmogorov-Obukhov  Theory 


The  above  discussion  has  been  relatively  general,  since  no 
specific  covariance  has  been  assvimed.  In  the  following  we  intro- 
duce the  von  Karman  spectrum  which  has  been  suggested  as  a pos- 
sible approximation  to  the  actual  wave  number  spectra  arising 
in  atmospheric  turbulence  [Tatarski,  1971] . This  spectrum  has  been 
verified  experimentally  in  certain  cases  but,  as  we  shall  see,  it  can  be 
necessary  to  use  correcting  terms  in  the  inner  scale  region.  The 
wave  number  spectrum  is  (nearly  homogeneous  and  isotropic  field 


assumed) 


-(f) 


r,  2 3 

o r 
n o 


* “TTi I 2 2 NTn/2  ' 3*  (2.8) 

n ^3/2  r(  ) (i  + 


The  corresponding  covariance  function  is 


V (r) 
n 


^2  (5-m)/2 

n 


r ( SCI  ^ 


(2.9) 
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In  these  expressions  is  again  the  variance  of  the  refractive 

index,  and  r may  be  considered  equal  to  the  outer  scale  L of 
o o 

turbulence.  The  exponent  m is  usually  called  the  refractive 
index  spectrum  slope  since  it  indicates  the  slope  of  the  spec- 
trum %<hen  plotted  on  a logarithmic  scale.  From  the  small 
argument  expansion  of  the  Bessel  function  K , v > 0,  it  can  be 

V 

verified  that  in  Bq.  (2.9) 

’’n  ■ 

Since  r(o)"a  it  is  seen  from  (2.8),  (2.9)  that  spectrum  slopes 
of  3 or  smaller  may  yield  negative  or  infinite  power  cp  (o)  1 
The  main  reason  for  this  is  that  the  correct  microstructure  of 
the  txirbulent  field  has  not  been  taken  into  account. 

If  the  von.Karman  spectrtim  is  used  a relation  can  be  found 

between  the  variamce  of  the  refractive  index,  and  the  structxire 

constant.  For  r « r it  is  fotmd  that 

o 


D_(r)  ~ 2 0. 


Iff)  ^ 


m— 3 _ _ 

r , 3 < m < 5 


iLx 

_/-l+m 


(2.10) 


For  this  to  be  valid  it  is  further  required  that  m < 5.  in  the 
special  case  of  m ■ 11/3,  as  predicted  by  the  Kolmogorov- 
theory  (see  Tatarski,  1971),  we  get 


2/3  ^ 2 
r C 
o n 


L(l) 

Ki) 


0,637  * 


(2.11) 


If  m ^ 5 it  can  be  shown  that 
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We  have  assumed  a locally  homogeneous  and  Isotropic  medium, 

which  is  invalid  when  r is  large  in  Eq.  (2.9),  or  when  k is  small 

in  Eq.  (2.8).  In  fact,  it  is  usually  possible  to  assume  that 

k r » 1 without  jeopardizing  the  validity  of  Eq.  (2.8), 
o 


....  2 , -m  3-m  ^\2j 

i..(k)  ^ o k r^  -77; 

n n o 3/2  _/m-3^ 


r(E=^) 


sin(^  (m-3)  ^ C ^ k”™,  3 < m < 5, 

4^2  V2  ) n 


(2.12a) 


In  particular,  for  m « 11/3  we  get 


♦ (k)  - 0.033  C ^ 
n . n 


(2.12b) 


This  is  the  form  in  whidh  the  spectrum  is  most  often  used  in 
turbulence  theory. 

As  noted  earlier,  if  k becomes  very  large,  it  is  necessary 
to  Include  the  effect  of  the  microstructure  of  the  turbulence 
(e.g.,  the  inner  scale  of  turbulence).  Tatarski  (1971)  has 
derived  the  expression  (for  m ■ 11/3) 


* (h) 
n 


0.033 


c 2 k-^^/3 


(2.13) 
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where 


k * 5.92/t  . 

m o 

The  approximate  stmcture  function  is 

r 

^ 2 . -4/3  2 

C -t  r 

n o 


D„(r) 


„ 2 2/3 

C r 
n 


r « K,, 


« r « 
o o 


y is  the  inner  scale  of  turbulence. 

'^o 

If  an  exponential  factor  like  the  one  in  Eq.  (2.13)  is  used 
for  other  values  of  m in  the  general  van  Karman  spectrum  then  the  restric- 
tion m > 3 is  no  longer  necessary.  However,  the  variance  of  the 
refractive  index  will  be  strongly  dependent  on  the  inner  scale, 
when  m < 3,  while  practically  independent  when  m > 3. 

2. 2. 3. 5 Theoretical  and  Empirical  Evaluation  of  the  Relevant 
Parameters 

In  Section  2.2.1  an  empirical  formula  for  the  mean  refrac- 
tive index  was  described.  It  relied  on  measurements  of  pressure, 

water  vapor  content,  and  temperature.  A similar  formula  can  be 

2 2 
found  for  the  variance  o or  the  structure  constant  C 

n n 

Tatarskl  (1971)  has  shown  that  the  following  expression  can  be 
2 

used  for  C in  the  Kolmogorov-Obukhov  theory. 


o 


(2.14) 


where 


c ^2.8, 


< 
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and 


I 


M 


79.10 


-6 


1 + 


15500q 

T 


Vdz 


V.  ' 


7800 


l,15500q 

T 


dz 


(2.15) 

In  the  above  expression, 

L Outer  scale  of  turbulence  (In  meters) 

o 

p Total  pressure  iri  m-bar 

T Temperature  in  degrees  Kelvin 

<I«0.62  e /P  Specific  humidity 
o 

e Water  vapor  pressure  in  m-bar 

o 

z Altitude  (in  meters) 

Y .0098^/Km  (adiabatic  temperature  gradient). 

d 

2 -2/3 

C then  has  the  dimension  of  m . Sirhis  (1971)  has  used  this 

2 

formula  to  determine  the  dependence  of  M on  altitude  and  humidity 

for  a typical  atmosphere.  The  result  is  shovm  in  Fig.  2.2.  Up 

2 

to  about  2 km  little  variation  with  height  is  found,  but  C can 

n 

change  two  orders  of  magnitude  when  the  humidity  changes.  Taking 

the  values  at  the  2 km  height  and  using  the  values  of  20,  50, 

2 

100,  150  m for  L yields  the  values  for  C listed  in  Table  2.1. 

o n 

TABLE  2.1 


Typical  Values  of  lO^^-C^^  (2  km)  Data  by  Sirkis,  1971) 


Lo/m 

Surface  Water  Vapor  Pressure  e^ 

0 

5 m-bar 

1 

15  m-bar 

30  m-bar 

20 

.68 

4.6 

38. 

106. 

50 

2.3 

15.5 

129. 

361. 

100 

5.8 

39.0 

325. 

910. 

150 

10.0 

66.9 

558. 

1,562. 
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Fig,  2.2  Structure  Parameter  as  Function  of  Altitude  with 
Sea-I«vel  Water  Vapor  Pressure  as  parauneter. 
(From  Sirkis,  1971) 
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In  Table  2 . 2 these  numbers  for  C are  converted  to  numbers  for 

2 " 
a using  Eg.  ( 2.11), 
n 

a ^ » 1.57  L C ^ 

n on 

2 2 

= 4.4  Lq 


TABLE  2 . 2 

2 13 

Typical  Values  Of  Refractive  Index  Variance  * 10 


L /m 
o 

Surface  Water  Vapor  Pressure 

0 

5 m-bar 

15  m-bar 

30  m-bar 

20 

.08 

.52 

4.4 

12.3 

50 

.5 

3.3 

27. 

77. 

100 

2.0 

13.2 

110. 

308. 

150 

4.5 

29.7 

247. 

593. 

To.  get  an  idea  of  typical  humidity  levels  we  list  below  the  water 
vapor  pressure  measured  by  radiosonde  at  Albany,  NY,  November  5, 
1976  (morning) . 

Altitude  (m)  Water  Vapor  Pressure  (m-^bar) 


86 

6.02 

145 

5.80 

256 

5.08 

733 

3.35 

927 

2.56 

1197 

2.84 

1297 

1.73 

1362 

3.48 

1456 

3.27 

1656 

3.81 

1995 

4.31 

2977 

2.73 
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Several  measurements  at  Albany  and  Buffalo  of  the  same  time-of> 
year  gave  very  similar  results  with  surface  water  vapor  in  the 
range  of  4-8  m-bar.  Lower  values  can  be  expected  in  colder  winter 
months,  while  considerably  higher  values  are  typical  in  the  summer. 

2 

M n __  _ _ . ^ ___  *•_  — _ ^ m .. 


Several  researchers  have  made  measurements  of  C 


Radar 


reflectivity  measurements  ^Kropfli  et.al.,  1968;  Hardy  and  Katz, 

2 

1969]  indicated  to  be  in  the  range 


Weak  scatter: 


Medium  scatter:  C 


20  • lO"^^ 

100  • 10“^^  m"^^^ 


_ 2 -«««  ,*-15  -2/3 

Strong  scatter:  6000  * 10  m 


The  two  first  numbers  match  fairly  well  with  the  results  in 
Table  2.1  for  50  m.  The  strong  scatter  returns  measured 


could  be  due  to  atmospheric  layering.  Most  <~75%)  of  the  reported 
measurements  lie  in  the  range  0-100  • lO”^^  . 


Broohner  (1970)  reported  on  some  measurements  of  the  optical 

2 

C performed  by  Goldstein,  et.al.,  (1965).  A typical  measurement 
” 2 

of  C vs  time  of  day  is  shown  in  Table  2.3.  Brookner  (1970) 

" 2 
then  uses  the  following  values  for  C as  typical  at  optical  wave- 
length 


Weak  scatter: 

c/ 

M .5 

10-^5 

n,"2/3 

in 

Med ium 

scatter: 

c = 

n 

~ 50 

10-^5 

-2/3 

m 

2 

-15 

-2/3 

Strong 

scatter: 

C 

n 

~ 500 

10 

m 

These  somewhat  arbitrarily  chosen  values  should  be  compared  with 
the  dry  air  column  in  Table  2.1. 


,**.*Wi*W«i»*^'**’* 
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I , L and  the  height  h are  all  in  meters.  As  an  example,  these 

O O j 

relations  yield  at  a 2 Ion  height,  j 


t ~ 13  mm 
o 

L ~ 90  m. 
o 

Superimposed  on  these  relationships  are  variations  with  time- 
of-day,  season,  weather,  etc.,  so  they  should  only  be  used  as  a 
general  guide. 

2 

The  variation  of  with  height  has  been  found  empirically 

to  be 

C ^ ^ ^ exp(-h/h_) 

n n no  o 

Pried  (1967)  used  b * 1/3,  h = 3200  m.  In  fact,  o ^ (h)  was 

o n 

found  to  be  of  the  form 

a^^(h)  - 6.7  • 10“^^  exp(-h/3200  m)  . 

These  results  are  for  optical  frequencies,  which  should  corres- 
pond to  dry  air  at  radio  frequencies.  A comparison  with  Pig.  2.2 
shows  that  the  expected  height  dependence  is  somewhat  weaker  than 
indicated  by  the  model  of  Pried.  On  the  other  hand  newer  data 
indicates  a much  stronger  height  dependence,  with  b - 2/3  and 
h^  320  m [see  Brookner,  1970,  1971],  Hence  the  issue  of  height 
dependence  is  not  quite  resolved,  although  it  appears  (see  Pig.  2.2} 
that  variations  in  the  first  2-3  km  can  be  ignored,  while  for- 
mulas of  the  form  above  can  be  applied  at  higher  altitudes.  Addi- 
tional data  on  both  the  optical  and  the  microwave  refractive  index 
variance  can  be  found  in  Gossard  (1977) . 

Prom  the  existing  data  we  can  conclude  the  following  typical 
values  (at  a height  of  1-2  km) 
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2 -2^3  -14 

C a.  10  m ' (t2ry  air)  to  10  fmaritime  air) 
n 

2 X4  12 

a 5 * 10  to  10  , depending  on  humidity; 

^ -14 

5 « 10  will  be  used  to 

represent  winter  conditions. 

L 70  m 
o 

I “*10  mm. 


2.2.4  Results  from  the  Theory  of  Turbulent  Scatter 

It  is  well  known  that  the  power  received  from  a small  scat- 
tering volume  can  be  written 


where 


P * P 

R T 


4ttR»j.^  4ttI^^ 


4tt 


(2.18) 


Transmitter  power 


Gain  of  transmitting  antenna 


Distance  front  transmitter  to  scatterer 


Scattering  cross-section  of  the  turbulent  volume 


R„  ■ Distance  from  scatterer  to  receiver 
R 


X » >iavelength 


Gain  of  receiving  antenna. 
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Eg.  (2.18)  also  serves  as  a definition  of  the  scattering  cross- 

section  (some  authors  define  a /4tt  as  the  scattering  cross-section) . 

s 

In  section  2.2.2  we  describe  some  general  relationship  derived 
for  atmospheric  layer  scattering.  Tatarshi  (1971)  has  derived 
the  cross-section  of  turbulent  scatter  with  a given  wave  number 
spectrum, 

a = 8ti  h $ (2k  sin  -f  ) sin  X dV,  i « X « L 

s n \ ^ / o o 


In  this  equation. 


(2.19) 


k * 2tt/X  is  the  wave  number. 


$^(k) 

n 

e 


is  the  locally  homogeneous  and  isotropic  wave 
spectrum  (e.g..  Eg.  2.8) 
is  the  scattering  angle 


sin  X accounts  for  loss  due  to  polarization  mismatch 
(usually  negligible) . 

X angle  between  incoming  electric  field  and  dir- 

ection of  propagation. 


dV 


is  the  infinitesimal  volume. 


Eg.  (2.19)  is  derived  under  the  assumption 


but  has  been  verified  at  much  higher  frequencies.  If  a spectnuD 

of  the  von  Karma'n  type  is  assumed  in  the  inertial  subrange  (i.e., 

K,  « X « L^) , it  is  found 
o o 


a 


s 


8 ^ k 


4-m 


^2  sin  "2  } * 


Polarization  losses  have  been  ignored  in  this  equation.  For 
m ■ 11/3  we  get  the  cross-section  prediction  by  the  Kolmogorov- 
Obukhov  theory,  while  for  m « 5 we  get  the  semi-en^irical  formula 
in  Eg.  (2.5).  The  data  for  Eg.  (2. 5)  are  mostly  at  frequencies  below 
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1 GHz.  Newer  experiments  at  higher  frequencies  indicate  that 
m = 11/3  provides  the  better  fit  to  the  data  (e.g. »Seehar8,1971) . At 
intermediate  frequencies  a mixture  of  the  different  effects  may 
be  described  by  an  equivalent  lii  in  the  range  3 < m ^ 5.  Combining 
the  above  expression  for  a^  with  Bq.  (2.18) 


R T T R 


r2,  ita,  t 


where 


_ 2 3— m , 2— m 

C=  CT  r k 
n o 


■(f) 


2M{^) 


(2.21)  i 


and  V is  the  total  scattering  volume.  Typical  values  of  the  at- 

2 . . 

mospheric  paxaaneters  and  r^  = were  described  in  the  previous 

section.  The  spectrum  slope  m is  11/3  for  turbulence  theory  and  1 

5 for  the  NBS  model.  Direct  measurements  of  m were  performed  by 
Eklund  and  Wicker ts  (1968)  by  comparing  the  received  power  at  1 
and  3 GHz  at  two  antennas  with  the  same  beeunwidth.  Average  values 
of  m in  the  range  of  11/3  to  5 were  measured,  and  a significant 
correlation  with  the  signal  level  was  observed.  This  is  consistent 
with  the  theory  that  both  turbulence  and  layer  reflection  are  im- 
portant at  these  frequecies  since  strong  layering  both  increase  m 
and  the  signal  level.  Eklund  and  Wicker ts  concluded  that  reflection 
is  dominating  at  1 GHz  suid  turbulence  at  3 GHz. 

Gjessing  and  McCormick  (1974)  also  report  measurement  of  m in 
the  range  of  2-6  and  give  an  empirical  foinnula  for  m (Gjessing  et.al., 
1969)  and  an  en^irical  distribution.  However,  the  validity  of  these 
results  is  apparently  limited  (Sherwood  and  Suyemoto, 1976)  particularly 
since  a strong  frequency  dependence  of  the  equivalent  measxired  m 

should  be  expected  from  the  discussion  eJsove. 
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2.3  Path  and  Antenna  Parameters 

In  this  section  we  list  the  parcuneters  that  are  required  to 
specify  the  troposcatter  path  completely,  and  briefly  comment  on 
the  significance  of  some  of  these  parameters.  Most  of  the  models 
developed  in  the  past  have  used  only  a sxibset  of  them,  or  have 
used  some  of  them  indirectly  through  other  auxiliary  parameters. 

A few  secondary  parameters  cheuracterizing  the  path  will  also  be 
discussed. 

2.3.1  Path  Geometry 

A typical  path  is  shown  in  Fig.  2.3  with  some  of  the  pcoram- 
eters  involved.  The  parameters  are 

• Frequency  f 

• Distance  d 

• Heights  h.p,hj^  of  transmitter  and  Receiver  antennas  above 
mean  Sea  Level  (MSL) . 

• Horizon  elevations  above  horizontal  trans- 

mitter and  receiver.  These  angles  can  depend  on  the 
atmospheric  conditions. 

• Antenna  boresight  elevations  above  horizontal 

transmitter  auid  receiver. 

Note  that  6qi-6«pG  and  9j(*-6j(G  then  the  elevation 
of  the  amtennas  adsove  the  horizon.  The  definition 
of  "boresight"  need  not  be  the  direction  of  maocimum 
gain,  but.  may  be  defined  from  mechamical  considera- 
tion (e.g.,  the  normal  to  a phased  array). 


• Azimuth  pointing  angles  cpip,cp|(  if  horizontal  diversity  is 
enmloyed,  or  beaun  swinging  experiment  performed. 


• Gain  patterns  Gt(P»cp)»  GR(9,cp)of  transmitting  and  re- 
ceiving antennas.  Determination  of  the  gain  patterns 
requires  a number  of  additional  parameters  - size  and 
type  of  antennas,  height  above  ground,  ground  profile 
and  reflectivity.  The  patterns  should  be  referenced 
with  respect  to  the  boresight  direction.  Both  anpli- 
tude  and  phase  of  the  voltage  gain  patterns  are  required 
when  the  cross  correlation  between  two  (diversity)  pathe 
is  to  be  calculated. 

• Polarization  of  transmitter  and  receiver 

• Bandwidth 

• Power 

2 . 

• Atmospheric  parameters  K,  9 , wind  velocity /attenuation, 

etc.  (See  section  2.2)  ^ 


The  performance  is  particularly  sensitive  to  the  antenna  eleva- 
tion angles  which  are  often  not  adequately  specified  in  actual  ex- 
periments . 

2.3.2  Antenna  Parameters 


In  most  theoretical  studies  it  is  assumed  that  the  antenna  pat- 
tern has  an  ideal  rectangular  shape.  This  ideal  pattern  is  extremely 
convenient  for  euialysis,  but  the  actual  pattern  is  substantially 
different.  An  equivalent  ideal  beam  shape  may  undoubtedly  be  postulated 
but  a general  relation  of  this  to  the  actual  beam  shape  is  not  well  defined 
Often  the  3 dB  beamwidth  is  arbitrarily  chosen  to  define  the  equivalent 
ideal  beam.  We  shall  later  (section  6)  estimate  how  good  such 
an  approximation  is  for  troposcatter  purposes.  For  sin^licity  the 
considerations  will  be  limited  to  the  parabolic  dish. 

The  gain  pattern  of  a parabolic  dish  can  easily  be  evaluated 
for  a uniform  illiimination  by  the  feedhorn,  as  well  as  for  s«ne 
specific  non  uniform  illvoaination  patterns. 
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Let 


i 


a - radius  of  the  circular  aperture 
A B Area  of  aperture 
X ® wavelength 
e ^ angle  off  bores ight 
r c distance  from  center  of  aperture 
b(r)  » illumination  of  aperture;  b(0)  * 1. 

The  voltage  gain  pattern  is  then 


g{e)  = cos  e “ J J dA 


■ 2ttx  . „ 
3—^  sine 


b(r) 


IM 


where  x is  the  distance  from  the  center  along  the  line  at  the 
intersection  of  the  plane  of  the  aperture  and  the  plane  determined 

by  the  incoming  ray  and  the  aperture  normal.  Hence » using  the  co- 
ordinates (r#  (p)  in  the  aperture# 

a 2tt  -Zrrr 

g('5)  = J rdr  J dcp  b(r)  e ^ 

•na^  o o 


coscp  sinP 


COsP 
2 

a o 


J dr  b(r)-  r sin  e)  . 


If  the  aperture  is  uniformly  illuminated  then  b(r)  « 1#  and 


g(6)  **  “2  COS  9 J 


If  we  tiry  a non  uniform  illximination  of  the  form 


0/(^  0 


b(r)  - (1  - (r/a)2)^  , 


than 


g(P) 


>M+1, 


g^(9)  ■ cos  6 2 (m+1)  1 
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J^^  sin  .) 


/2TTa  . - \ 

.ins) 


U+1 


The  half  power  beamwidths  are  determined  by  (ignoring  the  cos  9 
term) 


2TTa 

X 


sin 


M = 0 
M = % 
U = 1 


The  beamwidth  for  (i  = 1 is  approximately  23%  larger  than  for  uniform 
ill\imination,  euid  this  is  close  to  the  performance  that  is  obtained 
with  practical  dishes.  If  we  define  the  beamwidth  efficiency  of  the 
antenna  as 

h aurea  with  uniform  illumination  giving  actual  3 dB  width 
erriciency  - 


we  get: 


efficiency  = 


U » 0 

u * ^ 

M = 1 


It  can  also  be  useful  to  know  the  null-beamwidth  of  the  antenna. 

If  B ...is  this  beamwidth.  we  have 
null 


2TTa  / 

— V 


u = 0 

^ - Jj 

u “ 1 


The  nullwidth  with  non  uniform  illumination  is  much  wider  than  the 
3 dB  width.  Finally  the  level  of  the  first  sidelobe  can  be  found 
for  the  general  beam  pattern  above. 

-17.6  dB 

sidelobe  level  } -21.3  dB 

-24.6  dB 


u 


0 

*s 

1 
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For  the  remote  sidelobes  the  envelope  of  the  asymptotic  expansion 
of  the  Bessel  functions  can  be  used  as  a bound. 


^ 2 /Tia  . 

~ cos  6 -r-  sin  9 ) 

Tt  \ \ / 


G(e)  ® lg(9) ! ~ cos  e 

A nonuniform  illumination  is  often  used  on  purpose  to  achieve  the 
much  better  sidelobe  properties  at  a small  cost  of  aperture  effi- 
ciency. The  antenna  gain  6 can  be  evaluated  as  follows: 


4tt/G  * 


■i 


2tt  tt/2 

dcp  I de  sin  6 | 


Tl/2 


2 1+u  2 *^1+,. 

2tt  [ dP  cos  9 sin  8 4 ^(l+u)!  = — 


I 


(c  sin  e) 


2 (u+1) 


where 


Hence 


or 


c . » 1 

A 


_ - l+2u  r 2TTa-i 
® “ ^2  i IT 

A definition  of  gain  efficiency  analogous  to  beamwidth  efficiency 
results  in 

jloo%  U « 0 

gain  efficiency  ■ C 89%  U * *S 

I 75%  u - 1 
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The  Troposcatter  Model  with  Wide-Beam  Antennas;  Path 
Loss  Calculation 

The  path  geometry  of  Section  3 and  the  tropospheric  scatter 
model  of  Section  2.4  will  be  combined  to  find  the  performance  of 
some  typical « although  idealized,  troposcatter  links.  It  will 
be  assumed  in  the  following  that  the  earth  is  smooth  and  non- 
reflecting and  that  the  antenna  beams  are  pointed  just  above  the 
horizon. 

2.4.1  Theoretical  Development 

The  Integration  over  the  common  volume  will  assume  a spectrum 
of  the  von  Karman  type,  so  that  the  path  loss  is  given  by 
Eg.  (2.20). 


(2.22) 

where  C is  given  by  Eq.  ( 2.21).  The  integral  will  be  evaluate*'' 
by  integration  over  the  transmitter  and  receiver  elevation  angles 
e^,  Gp  and  the  distance  z perpendicular  to  the  great-circle  plane. 
The  integration  is  illustrated  in  Fig.  2.4  by  a view  in  the  great- 
circle  plane,  l.e.,  the  z-axis  is  perpendicular  to  the  paper. 


The  angles 

( 


®R  " 
°T  • 


oij^.  0^  are  defined 

- a/2% 

»_  + a/2R 

T e 


by 


where  is  the  effective  earth  radius  and  6^  range  from  the 
horizon  9^^)  and  up.  The  distance  d^  is 


d > 2 R sin  (d/2R  ), 

OB  B 


mk 
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and  the  scattering  angle  in  the  great-circle  plane  is 


The  infinitesimal  integration  volume  dV  is  projected  onto  the 

great-circle  plane  at  the  distances  from  the  transmitter 

or  OR 

and  receiver,  where 

R„  “ d sin  a„/sin  0 
OT  o ^ o 

R__  ■ d sin  a_/sin  0 
OR  o ^ o 

The  total  distances  from  the  point  of  integration  to  the  terminals 


/ 2 2 
R_  = y R„  + z 
T ^ OT 


2 . 2 


\ ~ ^OR  ^ 


The  scattering  angle  0 at  the  point  of  integration  can  be  deter- 
mined by 


2 


(2  .in  I) 


2 ^ (2  .In  -f)  * 2[  * ^Vl 

[(!«/)  (1*.^^)]*“ 


where 


z™  ■ z/R  _ 

T or 


*R  ■ ^/‘'OR  • 


It  is  now  assumed  that  the  strong  scattering  angle  dependence 

allows  us  to  ignore  points  where  z '>r>  R or  R _•  This  approxima- 

OR  OT 

tion  is  good  if  0 is  small  enough.  Hence  it  is  found 
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The  path  loss  in  Eq.  ( 2.22)  now  becomes 
tt/2  tt/2 

=2  J ■’“r  [ ('*ot'"’*or)  ®o(2 

°To  ®RO 


10 


-L/10 
P 


where 
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= - (i'  ¥) 


°To  ^'o  ®TG 


®Ro  “ ® 


RG 


and 


a = d/2R 
o e 


We  now  involve  the  small  angle  assumption. 


sin  e 


2 sin  “ 


R m + R « 
oT  OR 


d d. 
o 


L is  then  with  good  approximation  given  by 
P 


10 


-Lp/10 


1 r 
c — * 
^2  d J 


do,, 

a +e„ 


J '°T  °R> 


-m. 


o “TG 


n + 9 
'^O  RG 


d (ro-1)  (m- 


1-2)  ®RG  ®Tg) 


2-m 


C B 


(1 


d(m-l)  (in-2) 


’rg  " 


( 2.23a) 


In  this  equation  the  factor  C contains  terms  depending  on  frequency, 

spectrum  slope  m,  correlation  distance  r , and  the  refractive  index 

o 

variance.  From  Eq.  (2.21)  we  have 


2-34 


J 


2 3-m 

a r 
n o 


"(¥))• 


Assuming  for  simplicity  that  = ^tg~  ^ elevated  horizons) 
and  introducing  the  e3q)resei  on  for  C in  terms  of  the  knovm  pareun- 
eters  we  get 


3-lp/lO 


2 3-m  2-m  ^2-m 

4(m-l)(m-2)  '^n  o min 


( 2.23b) 


^e  approximations  made  are  good  when  d « and  m > 3.  Waterman 

(1958)  has  evaluated  the  integral  in  Eq.  ( 2.22)  exactly  for  m an 

even  integer  and  e„  = 0^^  * 0.  For  small  values  of  d/R  his 
’ TG  RG  e 

result  agrees  with  Eg.  ( 2.23). 

2 

Eg.  ( 2.23)  was  derived  assuming  constant  values  of  a and 

n 

r * L . The  empirical  results  of  Pried  (1967)  indicates  that 
©2  o 

a and  L depend  on  height, 
n o 


L - L (hA  ) 
o oo  o 


where  Fried  found  that 


3200  m. 


6.7  . 10 


L « 113  m 
oo 


The  discussion  in  Section  2.2. 3.5  Indicates  that  h may  be 

2 ® 
smaller  and  that  can  be  quite  different  from  the  values 

no 

indicated  above.  For  equal  height  of  the  terminals 
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h - R - (e  - P) 

min  RG  TG 


2^ 

/ 


(« 


-le  + e 

min  RG  TG 


)^)/ 


8P‘ 


min 


and  Lp  can  be  written  directly  in  terms  of  usually  specified 
parameter. 

2.4.2  Comparison  with  Other  Methods 

A number  of  empirically  based  methods  exist  to  predict  path 
loss.  Unfortunately  most  of  the  empirical  results  were  found  at 
frequencies  below  1 GHz<  where  layer  reflection  is  likely  to  be 
more  dominant  than  turbulent  scatter.  This  means  that  these 
methods  often  have  limited  applicability  at  higher  frequencies. 
One  of  the  most  widely  used  techniques  is  that  developed  by  NBS 
(Rice«  et.  al.,  1967).  At  distances  such  that  Pd  < 10  km  the 
predicted  path  loss  with  omnidirectional  antennas  is 

Lp  - 135.8+0.34  + 30  log  + 10  log  ^ + 30  log  6 88 

- -74.2  dB  + 0.34  • lo“^  Pd  + 30  log  f + 10  log  d + 30  log  6 

where  P is  the  scattering  angle  in  radians  and  all  other  units 
are  MKS.  This  has  the  same  form  as  Eq. ( 2.23)  if  m » 5.  This 
is  consistent  with  the  fact  that  the  NBS  method  is  mostly  based 
on  measurements  in  the  100-1000  MHz  range,  where  layer  reflec- 
tion is  expected  to  dominate  the  scatter  mechanism.  Hence  the 
2 

values  of  o and  r found  in  Section  2.2  for  turbulent  scatter 
n o 

are  not  directly  applicable, but  can  be  chosen  to  fit  the  above 
equation,  which  could  then  be  used  at  low  frequencies.  We  shall 

restrict  the  comparison  to  nearly  symmetric  links.  First,  let 

-3  2 

the  term  .34  * 10  Od  be  due  to  variations  of  with  height 

of  the  form 
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a 


2 


n 

Assvuning 


“ o ^ exp''-h/h  ] = exp{-d($+0  +6  )/8h  ). 
no  O T R o 

small  elevation  angles  it  is  then  seen  that 


10  1 
2.303  8h 

o 


or 

h = 1.6  km. 
o 

This  is  only  two  times  smaller  than  the  value  suggested  by  Fried  for 

optical  system, and  compares  favorably  with  other  data  (Brookner,  1971; 

2 2 

Gossard,  1977) . By  conmarison with  Eg.  ( 2.23)  a /r  is  found, 

no  o 

c o -2 

= 5.8  * 10  m 

h - o 


This  number  would  also  be  found  in  a humid  turbulent  atmosphere . 

At  higher  elevation  the  dependence  of  r^  on  h is  significant. 

This  is  reflected  in  the  NBS  method  by  the  use  of  a different 

formula  when  6d  > 10.  It  can  be  verified  from  the  above  numbers 

2 “2/3 

that  the  Fried  formula  for  9 r combined  with  the  rest  of 

n o 

Eg. (2.23b)  with  m » 5 deviates  by  less  than  1 dB  from  the  NBS 
method  in  the  altitude  range  300-3000m. 

We  can  now  compare  layer-ref ledtion  with  turbulent  scatter. 
Fig.  2. 5 shows  the  expected  path  loss  as  predicted  by  the  NBS  formulas; 
and  the  path  loss  is  predicted  by  the  turbulent  scatter  theory  with 
2 -14 

9^  * 5 • 10  (typical  for  dry  air) 

r ■ 70  m 
o 

m - 11/3. 

It  is  clear  seen  how  turbulent  scatter  is  dominant  at  higher 
freguencies  and  distances.  In  the  transition  region  (usually  in 
the  band  between  1 and  10  GHz)  both  layer  reflection  and  turbu- 
lence can  be  dominant,  depending  critically  on  actual  parameters 
and  therefore  also  both  on  location  and  season.  In  wet 
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air  the  turbulence  signal  can  be  20  dB  stronger . Several  methods  | 

other  than  that  of  NBS  are  discussed  by  Panter  (1972) . We  can  | 

mention  the  CCIR  method  (nearly  the  same  as  the  NBS  method),  Yeh's  I 

method  and  the  method  of  Collins  Radio.  | 

The  method  of  Yeh  predicts  the  path  loss  j 

Lp  » 30  log  -■  ^ + 20  log  + 573  d/R  - 0.2  (N  - 310)  + 57 

. >'  1 Mnz  X mi  e s j 

•i  t ^ 

■i 

where  is  the  surface  refractivity.  This  is  shown  in  Fig.  2.6  | 

for  a 5 GHz  link,  together  with  some  of  the  other  prediction  ] 

methods . | 

The  method  of  Collins  is  based  on  an  experimental  curve  at 
1 GHz.  An  excellent  approximation  to  this  curve  up  to  300  mi  ; j 

yields  1 

Lr.  “ 30  log  - ^ + 80  log  , + 34. 

P ^1  GHz  ^ 1 mi 

'Riis  curve  is  also  shown  in  Fig.  2.6  at  f = 5 GSlz.  The  received 
power  by  Collins  method  falls  off  much  more  rapidly  as  a function 
of  distance  than  does  any  of  the  other  methods. 

We  conclude  that  the  turbulent  scatter  theory  agrees  well 
with  other  methods  in  the  frequency  range  of  2-5  GHz,  and  that 
the  method  of  NBS  is  more  accurate  at  lower  frequencies.  This 

i 

means  that  lower  values  of  the  spectrum  slope,  and  higher  values 
of  the  refractivity  variance  should  be  used  at  lower  frequencies, 
but  the  model  with  a von  Kan^n  refractive  index  spectrum  may 
still  be  used. 
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Preliminary  CompariBon  With  anpirical  Results 

Numerous  measurement  of  tropes catter  loss  have  been  made  in 
the  past,  but  we  shall  only  chose  two  examples  for  conqparison 
here. 

The  first  example  is  a 210  Ion  European  path  at  12.3  GHz 
(Abel,  1971).  The  parameters  of  the  path  and  the  median  received 
field  are  listed  in  Table  2.4. 


TABLE  2.4 

Parameters  of  12.3  GHz  Link  (Abel,  1971) 


Transmitter : 

10  W 

Distance : 

210  km 

Frequency: 

12.3  GHz 

Total  Antenna  Gain,  G_  G_: 

T R 

80  dB 

Minimum  Scatter  Angle  (K  •=  4/3)  : 

11  m-rad 

Antenna  Diameter  (Both) : 

1 m 

Median  Field  Strength,  Winter: 

16  dB/-^ 
m 

Summer: 

22  dB/  ^ 
m 

The  antenna  beaunwidth  (33  m~rad)  is  large  enough  to  apply  the 
widebeam  results  of  Section  2.  4.1.  From  Table  2.4  the  measured 
path  loss  is  found: 

Path  Loss  (Winter) : 225  dB 

Path;  Loss  (Summer) : 219  dB. 
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The  "predicted  path  loss”  is  found  by  Inserting  the  above  para- 
meters in  Eg.  ( 2.23a),  using 


m = 11/3 

0^  = 5 • 10  (dry  air) 

and 

r = 70  m. 
o 

With  these  parameters  the  estimated  path  loss  is 

Path  Loss  (calculated)  « 220  dB, 

which  certainly  is  very  close  to  the  measured  loss.  If  a small 

coupling  loss  is  associated  with  the  path  the  theoretic  loss  will 

probably  be  closer  to  the  winter  measurements  where  the  humidity 

is  low.  The  stronger  signal  measured  in  the  summer  is  due  to  the 

2 

higher  humidity  (o^  larger,  see  Table  2.2). 

Now  consider  some  recent  data  on  a 4.8  GHz  link  (Sherwood 
and  Suyemoto,  1976) . The  parameters  and  measured  pathloss  on  this 
link  are  listed  in  Table  2.5. 

TABLE  2.5 

Parameters  of  4.78  GHz  Link  (Sherwood  and  Suyemoto,  1976) 

1000  W 

86  St.  mi. 

4.780  GHz 
78.2  dB 
8 ft. 

.75° 

-96  dBm 
-84  dBm 


Power,  P_; 

T 

Distance: 

Frequency: 

Total  Antenna  Gain,  G„  ■<-  G„: 

T R 

Antenna  Diameters  (Both) : 

Horizon  Elevation  (Both) : 

Median  Received  Signal  Level,  Winter: 
Median  Received  Signal  Level,  Summer: 


The  parameters  in  Table  2.5  yield  the  measured  path  loss 
or 


L = P„  + + G - RSL, 

\ p T T R 

Path  Loss  (winter) : 234  dB 

(svunmer)  : 225  dB  . 

This  is  now  compared  with  the  path  loss  predicted  from  Bq.  (2.23), 
using  again  the  parameters 

m = 11/3 

0^  = 5 *10  (dry  air) 

r = 70  m. 
o 

Inserting  the  parameter  above  and  those  of  Table  2.5  in  Eq.  (2.23a) 
results  in  the  estimated  path  loss; 

Path  Loss  (calculated)  » 222  dB. 

This  is  close  to  the  median  measured  in  the  summer.  Theoretically, 

with  the  assvimed  parameters,  the  results  should  be  closer  to  the 

winter  measurement.  However,  a small  coupling  loss  is  undoubtedly 

present  (5.2  dB  according  to  the  CCIR  formula)  and  the  uncertainty 

2 

of  the  atmospheric  data  (e.g., height  variation  of  9 and  r^) shows 
I no 

that  the  measurements  are  again  consistent  with  the  turbulent  scat- 
ter model.  In  Table  2.6  below  the  results  of  the  above  two  cam- 
par  isons  are  summarized. 

2.5  Spatial  Correlation  at  Receiving  Site 

The  derivations  in  this  section  follow  the  development  in  sec- 
tion 2.4.1  closely  and  we  will  therefore  refer  to  that  section  for 
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1 


Table  2.6 


Comparison  with  Empirical  Data  {Summary) 


Path  1 


Path  2 


Median  measxured  path  loss, 
winter  _ 


Median  meas.  path  loss, 
sximmer  


Path  loss  cai 
using  m=ll/3 
and  = 5 • 
coupling  loss 


; (Table  2.4)  (Table  2.5) 


225  dB 


219  dB 


234  dB 


222  dB 


222  dB 


intermediate  results.  Horizontal  and  vertical  correlation  will  be 
considered  separately.  The  case  of  a wide  beam  transmitter  is 
treated  first. 


2. 5 . 1 Horizontal  Correlation,  Wide-beam  Transmitter 


Suppose  a second  omnidirectional  receiver  is  situated  a dis- 
tance b from  the  main  receiver.  Assuming  b is  small  we  can  still 
use  the  integration  volume  in  Fig.  2 .4  for  both  received  power 
levels  as  well  as  for  calculating  the  correlation.  The  only  dif- 
ference from  the  integration  in  Section  2.4.1  is  that  the  inte- 
grand is  modified  by  the  phase  factor. 

Ids 

e 

where  z is  the  distance  from  the  great  circle  plane  in  Fig.  2.4. 
As  a reasonable  approximation  assume  that  z « bo  that 


2 ^ _2\h 


T- 


\ ^ • 
oR 
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The  correlation,  is  now  (compare  Eq.  (2.22) 


-J 


= c 


kbz 

oR 


ti  J 4I 

V 


Sin 


1\ 

2/ 


-m 


d^r 


The  integration  in  the  z direction  (perpendicular  to  the  paper  in 
Fig .2. 4 becomes 


where 


and 


. kb 


r 


“D 


dz  e 


oR 


(.  Sin  f)-  ^ 


2 _ 2 


o / 

= 2 „ 2 s 

OT  oR 


_ ^ 2.-nv/2 

2 sin  e (1+x  ) dx 


R ^ ^ 1-m 

o /.s  ON  _/l  m-l> 


2 2 

R m R « 
oT  oR 


(2sin-f)  B(f. 


e - ^ 2 R sin  ^ 
R „ o 2 

oR 


in 


(I) 


m-1/2 


Vl/2<»'- 


we  are  primarily  interested  in  finding  the  integral  of  t|  (b)  but 

Xl 

note  that  the  above  fonsula  can  also  be  used  to  find  higher  order 
moments  of  T\g(b).  If  we  write  0 ■ b * a*  where 

kR  e 

* * r?  ^ ' 

oR 
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we  use  the  well  known  Fourier  relation 


r 6(y) 


and  get 


J 


f r%  r 

a_ 

Ro  To 


1 -1 

J ^ — (sin  0 ) (2  sin  6/2) 


J 

TTC  r“  °Trs) 

kd  ^ 


-l-m 

- 6 . * 
ttC  min  _ 

kd  ra-1  , 


* 1 T ,K 

G J.  — I — f 2 \ 

^ k=0  Wo-*-  ®Ro^ 


In  this  equation. 


6 . * a_  + a„  = R + 0 + 6 

min  To  Ro  ' t TG  RG 


is  the  scattering  angle.  Let  us  define  the  function 


1 k 

Pn  (x)  * n r x , 0 < x < 1 

k«0 


the  following  bounds  on  are  useful 


(l  - X s e^(x)  s (l-x)“^. 


The  lower  bovind  is  often  an  excellent  approximation.  Thus  we  have 
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p 


CD 

J ri^(b)  db 
o 


0 . a 

ttC  inin  g / To  N 

^ “To-^  °Ro^ 


1 

2 


r,H(0) 


in-2 

m-l 


min 


/0  . ) 
To  min 

m-1^ 

2 / 


Ti  (0)  the  receiver  power  calculated  in  Section  2.4.1.  The  hor* 
H 

izontal  correlation  distance  is  defined  below. 


Definition  of  Horizontal  Correlation  Distance; 


Lr  = J db/r|jj(0). 


(2.24) 


With  this  definition  the  result  is 


n 3 , /9  . ') 

^ 1-2  m-l  \ To  min/ 

X “ m-l  „/l  in-lN 

^\2  ' ~) 


L„e  . 

H min  m' 


(2.25) 


Note  that  for  a symmetric  link  (e„  * 0__)  /0  . = 0.5  so  the 

TG  RG  TO  min 

right  hand  side  in  Eq.  (2.25)  is  then  independent  of  the  distance 
and  the  elevation  angles. 


2. 5 . 2 Vertical  Correlation 


We  now  determine  the  spatial  correlation  in  the  vertical 
direction.  Suppose  the  line  along  which  the  correla'tion  is 
bsken  is  tilted  a small  angle  & from  vertical.  6 determines 
how  far  the  cosibined  antenna  is  pointed  above  the  horizon. 
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In  this  case  the  phase  difference  between  the  two  received 
signals  is 


kb  sin  (6_-  6) 

Xx 


where  again  b is  the  distzuice  between  the  two  receivers.  The  cor- 
relation, is 


"jkb(6  -6)  -m 

- = Jj  J*  71— T (2  "i”  f) 


'*R 


Integration  along  the  z-axis  is  just  as  in  Section  .4.1,  so 


TT^(b)  = 


tt/2  ^Tt/2 

I dCL  dCL  e ^ ”(R  +R  „) 

2 J T J ^ R ^ oT  oR 


2 J T ' 
To  Ro 


0 « i-i-i 

sin  e (2  sin 

o 2 


S f ^ -^'^<V'^ro'*-®rg"^V 

d(in-l}  J \ 


+ 

To 


Here  is  again  the  constant 


_ _ _ 1 m-lN 

=2  " ®vl  ' ~)- 


Ti^(b}  can  be  written  in  terms  of  the  incomplete  Gamma  function  of 
a complex  argument. 


r(2-”'l“>(»To*“RoO  • 
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INriafcWi' 


L 


In  analogy  with  the  horizontal  diversity  the  correlation  distance 
should  be  defined  in  terms  of 


09 

J n (b)  db/ri  (0) 
J V V 


However,  in  this  case  it  is  important  that  the  aperture  is  aiming 
at  the  right  point  over  the  horizon, so  the  definition  must  include 

A 

an  optimization  of  the  angle  6, 

Definition  of  Vertical  Correlation  Distance; 


= max 
6 


r n (b)  db/n  (0) 


(2.26) 


It  can  be  verified  that  the  maximum  is  at  6 = 9__+  (cor- 

RG 

responding  to  physically  aiming  the  antenna  just  above  horizon) 
and  that 


maj^  n (b)db 
6 1-^ 


C 

2 In 

d(m-l)  K ^“"To  Ro^ 


so  that 


(2.27) 


Compared  to  the  horizontal  case  it  is  seen  then  is  laurger  than 

L by  a factor  of  2.6  - 3.4. 

H 

2.5.3  Comparison  of  Horizontal  and  Vertical  Correlation  Distance 
With  Results  by  Giessing  and  McCormick  (1974) 

G jessing  and  McCormick  (1974)  give  the  following  approximate 
expressions  for  the  correlation  distances  (also  assximing  a wide- 
beam  transmitter) : 
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0.44 

(4^”-  1)*“ 


0.44 
2^/”-  1 


These  formulas  clearly  assvime  that  the  horizon  emgles  eure  zero. 

A similar  method  by  Waterman  (1958)  yields  expresions  of  the  same 
form,  but  with  different  peurauneters . A comparison  of  these  ex- 
pressions with  ours  is  shown  in  Table  2.7. 


Table  2.7 

Comparison  of  Formulas  for  Correlation  Distance 
(Symmetric  Link.  Smooth  Earth) 


Spectrum  Slope  m 


(Eq.(  .25)) 

0.39 

.64 

0.65 

0.95 

1.18 

Gjessing  and  McCormick 

.57 

.65 

i .68 

.78 

.86 

(Eq.(  .27)) 

! 

1.00 

— 

1.67 

1 

t 

2.00 

3.00 

4.00 

> 

Gjessing  and  McCormick 

1.69 

2.11 

■ 2.33 

2.96 

3.59 

It  is  clearly  seen  that  there  is  only  little  difference  between 
the  results  here  and  those  of  Gjessing  and  McCormick  (1974).  The 
experimental  verification  in  their  paper  therefore  also  applies 
to  our  results.  However,  our  approach  is  more  accxirate  (asyn^)- 
totically  exact  integration) , more  general  (includes  elevated 
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horizon  angles) , euid  is  based  on  a definition  of  the  correlation 
distance  which  has  a direct  physical  meaning  (as  discussed  in 
Section  2.6). 


2. 5.4  Correlation  DistcUices  with  a Naurrow  Beam  TrcUismitter 

The  horizontal  correlation  fiuiction  is  extremely  wide  in 
this  case  since  the  receiver  only  sees  scattered  power  from  a thin 
vertical  column  in  the  horizon.  So, 


Tig(b)  ~ 1. 


(2.28) 


The  vertical  correlation  function  is  foxind  as  in  Section  2.5.2, 
but  integrating  only  in  a transmitter  beaun  with  the  solid  angle 

An 


OR  R 


R n R m 
oR  oT 


CfiQ  r “so*  p-n 


do^. 


Hence, 


l-m 

r,  (0)  .§42  ^ 

V d(X_  m-1 


and 


db  » X fi  , 

dou  min 


The  vertical  correlation  distance  is  then 
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X 


(m-l)  . 


( 2.29) 


OD 

L = J T)  (b)db/Tl  (0) 

V V V 


Pmin 


This  is  almost  the  same  as  for  a wide-becun  transmitter  (Eg.  (2. 24)). 
2. 6 Coupling  Loss  - Wid^eam  Transmitter 

Coupling  loss  may  be  defined  as  the  power  loss  incurred  by 
not  illuminating  the  scatterers  outside  the  common  volume.  Alter- 
natively it  can  be  defined  as  the  loss  due  to  the  decorrelation  - 
over  the  receiving  aperture.  These  definitions  are  equivalent 
when  the  transmitter  has  a wide-beam  antenna,  while  only  the  first 
is  directly  applicable  with  a narrowbeam  transmitting  antenna. 

Coupling  loss  is  a commonly  used  parameter  since  it  permits 
the  engineer  to  express  the  path  loss  in  terms  of  the  basic  path 
loss  with  omnidirectional  antennas,  the  antenna  gains,  and  the 
coupling  loss.  However,  for  an  actual  path,  such  a division  of 
the  path  loss  is  not  very  useful  since  it  requires  the  introduc- 
tion of  the  atmospheric  structure  outside  the  actual  antenna  becuns. 
Thus,  the  coupling  loss  will  depend  on  the  presence  and  strength 
of  inversion  layers  outside  the  common  volume!  Although  the  de- 
finition can  be  modified  to  be  based  on  a yearly  median,  say. 
this  would  imply  that  ^ort  term  variations  in  the  path  loss  has 
to  be  included  in  the  "basic  path  loss"  which  then  would  loose 
its  interpretation  as  the  loss  with  omnidirectional  antennas. 


The  logical  approach  is  to  calculate  the  path  loss  directly 
by  an  integration  over  the  common  volume,  using  the  formulas 
discussed  in  the  previous  sections.  This  integration  is  performed 
by  the  SIGHATRCan  TROPO-path  loss  computer  program. 
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1 

1 


For  historical  reasons  it  is  still  of  some  interest  to  calculate 
the  coupling  loss.  This  coupling  loss  will,  of  course,  only  have 
direct  physical  meaning  if  the  atmosphere  approaches  the  model 
atmosphere  assumed.  When  this  is  kept  in  mind  coupling  loss  can 
be  a convenient  analytical  tool. 

In  practice  a realistic  coupling  loss  is  not  well  defined 
in  the  1-10  GHz  range  where  layer  reflection  and  turbulent  scatter 
may  be  of  equal  influence,  depending  on  location,  season,  time-of- 
day,  etc. 

In  the  following  the  coupling  loss  for  a turbulent  atmosphere 
is  found . 


2.6.1  Coupling  Loss  for  an  Ideal  Narrow  Receiving  Beam 

Consider  a narrow  receiving  beam  with  unit  gain  in  a small 

solid  angle  and  zero  gain  outside.  It  is  assumed  that  all  of 

this  angle  is  pointing  over  the  horizon.  If  L is  the  path  loss 

pn 

with  this  narrow  beam,  and  Lp  is  the  widebeam  path  loss  formed  in 

Section  2. 4.1, the  coupling  loss  l (in  dB)  is 

c 

L = L - L . 
c p pn 


Lpn  is  now  calculated.  The  common  volume  is  a section  of  a cone, 
but  we  can  still  reference  to  Pig.  2.4  for  the  integration. 
is 


10"^Pn 


10 


(2  sin 


m 

) d 


tt/2 

» can  J dd^ 

°To 


1 1 
«oT  ®o 


(2 


sin 
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Aaaiiming  that  d « R , thla  baconas 

6 


=.  — SAG 


(<^e  °To)^'" 


d sin  m-l 

/A  + 0 +0  >1"® 

\R  ®RG  ®TG' 


C£SL 
d 


‘^(2R  ®Rg) 


m-l 


From  this,  and  Eq.  ( .23a)  the  coupling  loss  is 


Lc  = 10  log 


/'A 

VR 

e 


8rg  «tg)  ( 25"  «rg)  ®(i' 


AO 


m-2 


( 2 .30) 


This  expression  is  only  valid  asymptotically,  e.g.,  for  small, 
solid  angles  AO*  If  the  beam  is  rectangular  with  horizontal  and 
vertical  bean»/idth  and  P^,  respectively,  then 


■ «h  • Sv 

For  a circular  antenna  with  bearawidth  r,  ao  is 
AO  “ TT  (r/2)^. 

For  a linh  where  6^^^  the  coupling  loss  depends  directly  on 

e^in/A0»  where  is  the  angular  distance  (minimum  scattering 

angle) . 


2.6.2  Coupling  Loss  of  a Uniformly  Illuminated  Receiving  Aperture 

The  coupling  loss  can  also  be  determined  from  the  correlation 
distances  calculated  in  Section  2.5.  Suppose  E (r)  is  the  field 

O 
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received  at  a point  £ in  the  receiving  aperture  A . Assuming 
that  E(r)  is  locally  homogeneous,  let  p (£  - £' ) be  the  correla- 
tion  function, 

p(r  - r')  = e{  E^(r)  * E^(r*)}  • 


The  received  power  is  then 

■ jJ  <><£-£•> 


Perfect  correlation  would  have  yielded  the  maximum  power. 


R,  max 


2 

P . 


Assume  for  simplicity  that  p(o)  = 1.  This  entails  no  loss  of 
generality.  The  coupling  loss  is  therefore 

lO-VlO  . _l_  IJJJ  ar  dr'  p(r  - r') 

*R 


Exact  calculation  of  is  not  possible  unless  the  complete 
correlation  function  p(r  - x')  is  known,  but  simple  approxima- 
tions often  yield  nearly  the  correct  result,  consider  first  a 

rectangular  aperture  of  dimensions  a , a , so  that 

H V 


The  coupling  loss  can  then  be  written  in  the  form 
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10 


-Lo'lO 


a a 
H V 


I 


dy 


a. 


:r)  ^ 


-a 


V 


where  x and  y here  are  the  horizontal  and  vertical  aperture 
coordinates,  respectively.  The  following  approximation  is  often 
adequate: 

p(r  - £*)  = p(x  - x',  y - y*) 

- px  py  (y  “ y’)  • 


This  yields 


If  only  the  asymptotic  coupling  loss  is  considered,  i.e.,  a , a 

H V 

large,  then 

*H  ®V 


where  and  are  the  correlation  distances  calculated  in 
Section  2.5.  Using  Eqs.  ( 2.25),  ( 2.27)  we  get  the  approximate 
coupling  loss  for  rectangular  aperture: 


L -10  log 
c ^ 


2 

(°To^®min) 

a a 

H V 

®TG  ®RG/  , _ 

e (m-1)  B 

(i 

m-l\ 

2 ) 

( 2.32) 
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If  the  separability  of  the  horizontal  and  vertical  coor- 
dinates is  not  assumed,  the  following  more  accurate  result  is 
obtained  by  integration  over  the  common  volume. 


( 2.33) 

This  is  the  same  as  Eq.  ( 2,30)  if  the  solid  beam  angle  of  the 
equivalent  ideal  beam  is 


( 2.34) 


where  A is  the  aperture  area.  I^is  relation  assumes  equal  gain 
for  both  the  aperture  and  the  ideal  beam. 

A comparison  of  Bqs.  ( 2-3Z  and  ( 2.33)  shows  that  the 

separability  assumed  in  Eq.  ( 2.32)  leads  to  slightly  different 

values  for  L . However,  for  a link  where  ■ 0^  the  difference 
C R6  TG 

in  coupling  loss  is  only 

‘S  - 1“  log  ^ Vi'**') 

which  is  in  the  range  -4.1  dB  to  0 dB. 


2-57 


Let  us  next  consider  a circular  aperture  with  diameter  D. 

For  the  purpose  of  determining  the  asyn^totlc  coupling  loss  It 
Is  easy  to  see  that  only  the  area  of  the  aperture  Is  significant 
(this  assumes  that  p (£)  decreases  sufficiently  In  all  directions). 
Then  coupling  loss  for  circular  aperture; 


L = 10  log 
c ’ 


ttD 


4x  e 


— + e__  + + s— 1 


®TG  "'RGA2R 


1 m-l> 

e 

m-: 


( 2.35) 


The  beamwldth  of  the  equivalent  Ideal  circular  beam  Is  simply 


= - = 1. 


p = - = 1.27 

TT  D D 


Fbr  a uniformly  Illuminated  aperture  this  Is  found  to  correspond 
to  the  4.8  dB  bearorfldth,  using  the  results  of  Section  2.3.2. 


2.6.3  Comparison  With  Other  Results 

For  a link  with  zero  horizon  elevation  angles,  G jessing  and 
McCormlcK  (1974)  suggest  the  formula 


^2  2 

* , (L  ) . 5 (4^^”  -I)*®. 

Ad^  = 


^ (L^)  not  in  dB^.  Waterman  (1958)  derives  a similar  expression: 


■J 


(5!/“  .1)  (s^/ro  . 


where  A Is  the  beamwldth.  From  the  discussion  of  Section  2.3.2 
and  the  results  of  the  previous  section,  a reasonable  expression 
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for  m is 


R 


1 

D • 


The  above  mentioned  results  are  compared  with  ours  in  Table  2.8. 


TABLE  2.8 

Comparison  of  Formulas  for  Coupling  Loss 


Spectrum  Slope  m: 

3 

11/3 

4 

5 

6 

X" (L^) 

Eq.  ( -35) 

1 

. 0.47 

0.39 

0.22 

0.15 

G jessing  and 
McCormick 

1 

0.71 

0.61 

0.42 

0.31 

Ae  . ^ 
mm 

Waterman 



0.98 

0.65 

0.55 

D 

0.26 

The  essential  difference  between  these  expressions  is  a 
slightly  stronger  dependence  on  m in  our  result.  Also  the  ad  hoc 
methods  of  G jessing  and  McCormick,  and  Waterman,  seem  to  over- 
estimate the  coupling  loss  by  up  to  3 dB. 

The  coupling  loss  formula  in  Eq.  ( 2.35)  can  also  be  compared 
with  experimental  data.  Consider  the  path  used  by  Crawford,  et.al., 
(1953) . nie  parameters  for  this  path  are  listed  below: 


Frequency  f: 

Transmitting  Dish: 

Receiving  Dish,  D: 

Angular  Distance 

Coupling  Loss  (relative  to  8 ft. 
Dish)  Measured  as: 
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4.1  GHz 
10  ft 
60  ft 

0.9°  (k  - 4/3) 


11  dB. 


In  setting  these  parameters  into  Eq.  ( 2.35),  using  m « 11/3, 
yields: 

Calculated  asymptotic  coupling  loss:  7.6  dB. 

The  difference  is  due  to  the  fact  that  the  asymptote  has  not 
yet  been  reached  (see  Fig.  2,7  for  the  typical  behavior  of  the 
coupling  loss) , and  that  the  transmitting  beam  is  not  exactly 
omnidirectional.  In  Section  2.7  we  shall  attempt  to  find  how 
the  actual  coupling  loss  approaches  the  asyn^totic  results. 

2.7  Coupling  Loss  for  Narrcw-Beam  Antennas  at  Both  Sites 

In  the  previous  section  a widebeam  transmitter  was  assumed. 
On  most  practical  paths  the  coupling  loss  is  associated  with 
narrowing  of  both  transmitting  and  receiving  beams.  In  this 
section  the  case  of  extremely  narrow  beams  of  comparable  size 
are  considered. 


2.7.1  Coupling  Loss  for  Ideal  Narrow  Beams 

The  loss  is  dependent  on  the  relative  horizontal  beamwidths 
of  the  two  beams.  Let  the  beamwidths  of  two  rectangular  beams 
be 

e„  ,,  = vertical  beamwidth  at  transmitter, 

T,  V 

5^  g “ horizontal  beamwidth  at  transmitter, 

B * vertical  beamwidth  at  receiver, 

R/  V 

P.,  „ * horizontal  beamwidth  at  receiver. 


Assume  that  the  receiving  beam  is  the  narrowest  horizontally. 


^R,H  ^oR  ^T,H  ^oT 
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If  this  is  not  the  case,  the  subscripts  T and  R should  be  inter- 
changed in  the  remaining  expressions  in  this  section. 

Since  the  integration  (see  Fig.  2.4)  is  now  over  a very 
small  volume,  it  is  simple  to  write  down  the  path  loss  Lp.  Using 
Eg.  ( 2.22)  we  get 


j-Lp/10 


- 1 „ -m  ‘ ^R^R.V  ' ^^T.V 

R 2 tnin  * 

R T 


d '’r,V  ^R,H  *^T,V  ®min 


-1 

°Ro 


Combining  this  with  Eg.  ( 2.23a)  yields  the  coupling  loss 


10  log 


* °RO 


e . 

min 


®R,V  ®R,H  ^T,V 


(m-l) (m-2) 


( 2.36) 


where 


®min  * R ®TG  ®RG 

e 


°Ro  “ 2R 


Loosely  speaking,  the  coupling  loss  depends  on  (o/b)  for  two 

2 

egual  narrow  beams,  and  on  (e/B)  only  one  narrow  beam 

(Eg.  2.30).  Viewed  as  a function  of  freguency,  the  coupling 
loss  has  a slope  of  9 dB/octave  for  two  narrow  beams  and  6 dB/ 
octave  for  one  narrow  beam. 

Following  the  general  principles  of  the  preceding  section, 
the  coupling  loss  can  also  be  expressed  in  the  form 
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■-  u/—  J-N 

V 2 ®V2*  2 ' 

10  log  ‘^Ro  (m-l)  (Tn-2)  (2.37a) 


a™  r.  / < R « / R (2.37b) 

T,H  R,H  oT  OR  ' 

where 

® Area  of  transmitting  aperture 
a ® Vertical  dimension  of  receiving  aperture 

R ^ V 

H^R  H * Horizontal  dimensions  of  transmitting  and 
* ' receiving  apertures,  respectively. 

When  inequality  in  Eq.  ( 2.37b)  is  reversed,  the  role  of  the 
transmitter  and  receiver  should  also  be  reversed  in  Eq.  ( 2.37a). 


2.7,2  Coupling  Loss  With  Intermediate  Values  of  Vertical  Beam- 
widths  . 

So  far  we  have  only  considered  the  asymptotic  case  of  very 
narrow  beams.  Here  this  assumption  is  relaxed  in  the  vertical 
direction,  while  the  horizontal  beamwidth  will  be  assumed  either 
very  large  or  very  small.  The  ideal  beamshape  will  be  used, 
i.e.,  the  beam  pattern  is 


G(e) 


fi  I el  < b/2 


0 1 el  > p/2 


where  B is  the  beamwidth. 

If  the  beam  is  wide  horizontally,  then  the  path  loss  is 


(from  Section  2.4.1) 

^„/l  m-l\  ®RO  ^RV 


,-Lp/lO 


2 > 


°to'^‘*tv 


‘W 
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where 


^(2*  2 ) 2-in  p/frv  fsv  N 

d(m-l)(in-2)  min  ' 


(Xj^,  x^)  * (1  “ (1+Xj^)^”™  -d+x^)^  ”*  +(1+Xj^+X2)^ 


The  coupling  loss  is  simply 


- -10  log  '■(''T/Un'  "R/Vin> 


( 2.38) 


For  wide  vertical  beams  this  reduces  to  Bq.  ( 2.23a) » while  for 
narrow  vertical  beams  we  can  use 


F(Xj^,  Xj)  =-  (m-1)  (m-2)  Xj^  x^. 


If  we  Instead  assiune  that  one  of  the  beams,  the  receiving  one  say, 
is  narrow  horizontally,  then  we  have 

^ “ro'^^rv  °To''''^TV 


°Ro‘^^RV 


dfin^Lfl  ^°R  a«  (“r'^°To)  "(°r‘^°To'^^Tv) 
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CP 


ffl—  ^ ^ /_f22 Nj 


\l-m  /°To'*'^TV 


Nl-m  /”ToTV  \ 

(®min  ®Tv/  (e  . +R.„J 


min  TV 


. / ^ . \l-m  / °^To^  ^TV  \ 

\®min‘‘‘^Tv‘*'''Rv)  ^m-1  Ve  . +«-,,+ Bo„) 

min  TV  RV 


The  function  6 ■.  (x)  was  defined  in  Section  2.4,  EUid  is  a special 

m—1 

case  of  the  hypergcometric  function. 


9 (x)  = F(v,  1;  V + 1:  x) 

V 


The  above  expression  for  the  path  loss  may  be  easier  to  interpret 
if  written  as  an  infinite  sum. 


10 


-Lp/10 


eg 


RH 


k«0 


d (m-1)  (m-l+k) 


Ic  „l-lc-m  \l-b-m 

°To  ®min  °To  V®min  ^Rv) 


The  coupling  loss  is  then 


CD 

L « -10  log  ^7 r * Z — 77.  , . 

C ^ . „/l  m-l+k 


^°To  \ 
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where 


.1-m-k  / 

e . k 

. mxnx  . 

^1+X 

V RVy 

' " V 

TV  Oto  ' ' 

^ TV^' 

6 . k 

1-m- 

( 

( i+x_  +X  ^ 

> 

\ TV 

“to  ^ 

\ TV  RV/ 

/ 

l~m-k 


( 


2.39) 


^TV  “ ^TV^®min 
\v  “ ®RV^®min  * 


This  is  shown  in  Fig.  2.7  for  a particular  link,  together  with 
the  asymptote  , Eq.  ( 2,36)  and  the  curve  for  a wide  horizontal 
beam,  Eq.  ( 2.38).  Combining  the  curves  for  Eqs.  ( 2.38)  and 
( 2.39)  yields  an  excellent  approximation  to  the  coupling  loss 
when  this  is  less  than  5-7  dB.  When  the  coupling  loss  is  higher 
Eq.  ( 2.39)  is  sufficiently  accurate  by  itself. 

Although  Eq.  ( 2.39)  is  not  difficult  to  evaluate  numerically 
it  may  be  worthwhile  to  note  that  the  following  approximation  may 
be  used  at  values  of  the  coupling  loss  larger  than  apprracimately 
10  dB. 


~ —10  log 


RH 


®Ro®b'  2 ) 


where  P(*  , •),  and  were  defined  above.  Numerical  calcula- 

tions indicate  that  the  error  is  less  than  1 dB. 
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The  derivation  of  Eq.  (2.39)  assximes  that  one  beean  is 
narrow  horizontally  in  comparison  to  the  other  beeun.  It  has 
been  verified  niimerically  for  symmetric  links  where  the  two 
beauns  are  equals.  In  the  importamt  case  of  nonsymmetric  links 
with  equal  horizontal  beamwidths  (assumed  narrow  as  in  Eq. (2.33)) 
the  following  approximation  and  lower  bound  is  easily  found:* 


- - 10  log 


2(m-2) 


H 


mB 


/I  m-l\  8 
\2'  2 / 


min 


^m-1  ^v) 


+ 10 


log  ( 


1 + 


min 


+ 8 - 8 
RG  TG 


I) 


(2.40) 


where 

Vi  - 1 -(l  * - (i  + 

and 

hj^  - h^  = difference  in  elevation  of  receiver  and 
transmitter 

difference  in  elevation  angles  of  receiving 
and  transmitting  auitennas. 

This  formula  reduces  to  Eq.  (2.36)  in  the  asymptotic  case  of 
very  narrow  beauns. 

2.8  Conclusions 

The  preceding  analysis  leads  us  to  the  following  conclusions: 
* The  concept  of  coupling  loss  should  not  be  used  as  a 
universal  answer  to  the  question  of  antenna  gain  loss 

* The  Beta  ftinction  *1*68  for  m * 11/3 
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1 

■ ^ 


on  a troposcatter  channel.  The  total  path  loss  is  a 
more  fundamental  parameter  since  it  only  depends  on 
the  atmosphere  inside  the  common  volume. 

[ 

• Coupling  loss  should  only  be  used  in  connection  with 
a specific  model  of  the  atmosphere.  Hence,  it  makes 
little  sense  to  compare  the  coupling  loss  for  different 
atmospheric  models.  In  particular,  the  coupling  loss 
we  have  calculated  for  m « II/3  cannot  be  compared  to 
that  of  NBS  {m  - 5)  nor  to  that  of  CCIR. 

• The  pathloss  can  be  calculated  directly  from  the  atmos- 
pheric pareuneters  of 

effective  earth  radius 
variance  of  refractive  index 
outer  scale  of  turbulence. 

This  provides  a direct  physical  background  for  finding 
the  long  term  distribution  of  the  received  power  level, 
and  opens  the  way  for  more  accurate  prediction  of  the 
performance  of  specific  links. 

• The  turbulent  scatter  theory  is  accurate  at  higher 
frequencies  (>  1 GHz)  while  layer  reflection  can  be 
dominant  at  lower  frecpiencies.  The  simultaneous 
existence  of  the  two  scatter  modes  can  complicate 
accnirate  link  prediction  at  all  times,  but  turbulent 
scatter  represents  worst  case  and  should  therefore 
always  be  used  for  link  design. 
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• For  best  results,  the  pathloss  should  be  calculated 

by  numerical  integration  over  the  common  volume,  using 
the  actual  antenna  patterns. 

• Highly  accurate  formulas  have  been  derived  for  the 
pathloss  for  the  turbulent  scatter  model.  The  pathloss 
is  expressed  in  terms  of  the  basic  pathloss  (Eq.  2.23a) 
with  omnidirectional  antennas,  and  the  coupling  loss 
(Eqs.  2.38,2.39,2.40) . The  accuracy  of  the  coupling  loss 
will  be  verified  below,  along  with  a summary  of  the 
relevant  formulas. 

• The  computer  program  and  the  formulas  derived  yield 
results  for  the  pathloss  which  are  very  close  to  those 
predicted  by  the  NBS  method  when  a spectrum  slope 

m a:  5 is  used. 

The  relevant  formulas  and  parameters  needed  are  listed  below: 
The  Scatter  Model: 

Von  Karman  Spectrum  with  the  spectrum  slope  parameter  m, 

3 < m. 

The  resulting  scattering  cross-section  is  described  in  Section 

2.2.4. 

The  Basic  Pathloss  (Omnidirectional  antennas)  : 


* 6.7*10”  ^exp(-h/3.2  km),  h is  the  height  of  the 
scattering  point.  This  is  only  to  be  considered 

- as  a typical  exzunple  since  considerzd>le  variation 
with  climate  and  location  exists  (Gossard, 1977) 

r^  * correlation  distance  of  the  refractive  index  fluctua- 
tions. For  turbulence  this  may  be  equated  to  the 
outer  scale  of  turbulence  Lq.  With  the  results  of 
Pried  (1967)  this  yields 

» L * 2(h/lm)^. 
o 

h “ 2tt/x  “ 2TTf/c 

* wavenumber  of  the  link  frequency. 

d • path  distance. 

R * effective  earth  radius 

e 

* horizon  elevation  relative  to  horizontal  at  receiver. 

•=  horizon  elevation  relative  to  horizontal  at  transmitter. 


The  Coupling  Loss 

Several  formulas  have  been  derived  for  the  coupling 
loss.  One  represents  the  coupling  loss,  ® wide  hori- 


zontal beam,  and  the  others  the  coupling  loss, 
horizontal  beam. 


with  a narrow 


The  formulas  assume  that  the  antenna  beams  have  been  approxi- 
mated by  ideal  rectangular  beams.  The  validity  of  this  approximation 
was  carefully  analyzed  for  the  asymptotic  case  (large  coupling  loss). 


The  results  (Eqs. 

2.38  -2.40) 

are 

fl  + 

\2-m 

- (i* 

®RV  N 

2-m 

1 

■ «n.in 

I 

®min^ 

1 

. f % 

+ ^ 

* «TV 

2-m 

®min 

) 

2 

-70 

where 


9 . ^ minimum  scatter  angle 

min 

= angular  distance. 

= vertical  beamwidth  at  the  receiver 
X/D  in  radians. 

0^  = vertical  beamwidth  at  the  transmitter 

=“  X/D^  in  radians . 


For  L_,„  two  results  were  obtained.  If  one  beam  is  narrower  than 

v-Nxi 

then 


. 1 m-2  ^RH 

^CNH  - - 10  109  ,i  5s^,  • "i • S, 


B 


mm 


where 


B(x,y)  is  the  Beta  function. 


00  1 / °To  \^  , /,  . ®RV  Nl-m-k 


k*o 


min 


^ 1-m-k 

. + -TV')  A + ^ 'I 

^o  ®min 


e„ 


0™.  k 

.(1*^)  (1.  - 


^ l-m-Kj 
min  J 


and  where  we  have,  in  addition  to  the  parameters  defined  above. 


B horizontal  beamwidth  at  the  receiver,  assumed 

narrower  than  at  the  transmitter, 

d 

<ipQ  * 2^  at  equal  height) 

e 
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r 


If  the  horizontal  beamwidth  of  the  two  antennas  are  identi- 
cal then  the  loss  L is  bounded  by,  and  approxi- 

RH  Tn  B CNxi 

mately  equal  to  for  very  narrow  beams. 


L > - 10  log 

caiH  - ^ 


2 (id-2) 


H 


„/l  m-l\  e . 

."Y";  min 


min 


min 


_TJ 

e 


min 


+ 10 


log^ 


1 + 


To  Ro 


I/e  . ) 

min/ 


(2.40a) 


where 


a - a 

RO  To 


2(V  h^) 


d 


+ (6^^  - 

RG  TG 


and  h^-  is  the  height  difference  of  the  antennas. 

To  test  the  accuracy  of  these  analytical  formulas  they  have  been 
compared  to  the  results  of  a direct  integration  in  Fig.  2.8.  nie 
stipled  curve  in  Fig.  2.8  has  been  found  by  combining  from 


Eq.2.39  and  by  the  expression 

CWxi 

& 7 log  [10*'=™'’  7 10''=™^'' 


] . 


The  form  of  this  equation  was  chosen  arbitrarily  to  get  a smooth 
transition  between  the  results  for  a narrow  and  a wide  horizontal 
beam.  The  number  7 was  chosen  to  fit  the  results  of  the  numerical 
integration.  Perfect  agreement  is  then  found  for  coupling  losses 
larger  than  5 dB.  For  smaller  coupling  losses  the  numerical 
integration  routine  becomes  inaccurate  due  to  the  large  common 
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volume  so  it  appears  likely  that  the  expression  for  L is  ac- 


curate even  in  this  range. 


The  formula  (2.40)  for  equal  size  antennas  can  be  adapted 
to  work  over  the  entire  range  of  beamwidths  by  using  the  fol- 
lowing expression 

P_..  P. 


...  _ /TV  RV  N 

Lc  * - 10  VAe—  ' — ) 


e . 

mm  mm 


10 


log( 


1 + 


« . 1 m— 1.  . 

2 2 mm  A 

e / 

H 


2(m-2) 

+ 10  log(l  + |a„  - 1/6  . ) 

To  Ro  mm  / 


(2.41) 


This  formula  is  asymptotically  correct  md  is  an  excellent  ap- 
proximation for  intermediate  size  antennas,  as  evidenced  by 
the  points  calculated  in  Fig.  2.8. 


1 


Fig.  2.8  Comparison  of  Derived  Coupling  Loss  Formulas  with  Numerical 
Integration  Results  for  A Symmetric  Link  (100  mi.,  1 GHz, 
smooth  earth , turbulent  scatter) 
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SECTION  3 


ANGLE  DIVERSITY  PERFORMANCE 

In  this  section,  the  troposcatter  model  for  performance  pre- 
diction is  used  to  design  angle  diversity  systems.  A qualitative 
approach  is  followed  to  provide  general  guidelines  for  angle 
diversity  implementation  and  three  specific  examples  are  used  to 
provide  a quantitative  base.  The  three  excunples  are  the  C-band 
RADC  troposcatter  test  system,  an  operational  DCS  L-band  tropo- 
scatter system,  and  an  operational  NATO  L-band  troposcatter  sys- 
tem. The  parcuneters  of  these  links  are  summarized  in  the  follow- 
ing t2ible.  This  section  also  presents  the  long  term  variability 
model  to  be  used  with  the  troposcatter  model  for  prediction  of 
link  availability.  The  question  of  a long  term  decorrelation 
advantage  of  angle  diversity  vs.  space  diversity  is  also  discussed. 

3.1  Angle  Diversity  Design 

: i 

The  short  term  (~  seconds)  fading  characteristics  of  tropo-  ' I 

scatter  signal  reception  are  mitigated  by  the  provision  of  redun-  i 

dant,  i.e.,  diversity,  channels  and  an  associated  combining  scheme  t 

i 

to  utilize  the  stronger  received  signals  and  eliminate  the  weaker 
ones.  One  method  of  realizing  additional  diversity  channels  is 
to  use  multiple  feedhorns  at  the  focal  point  of  a parabolic  re- 
flector to  realize  multiple  angle  of  arrival  signals.  This  method 
is  normally  only  employed  at  the  receiver  since  the  production  of 
multiple  transmit  beams  requires  either  additional  power  amplifiers 
or  a reduction  in  power  with  power  splitters.  Thus,  conventional 
angle  diversity  has  one  transmit  beam  and  two  or  more  angle  of  ar- 
rival receive  beeuns. 
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Table  3.1 
Path  Pareuneters 


Name 


Path  Distance (mi. ) 


Transmitter  Height ( f t . ) 


Receiver  Height (ft.) 


Transmitter  Antenna (ft. ) 


Receiver  Antenna (ft.) 


Transmitter 

Boresight/Horizontal  (°) 
Boresight/Horizon  4 (O) 


Receiver 

Boresight/Horizontal  (°) 
Boresight/Horizon  i (°) 


Squint  Angle  (®) 


Frequency  (GHz) 


rADC  S.Tepesi 
Teat  Link!  Yamanlar 


Kristiansand 


168.3 


340. 


460. 


28. 


1.27 


0.77 


0.27 


0.54 


.5 


2. 


2893. 


3176. 


60. 


60. 


0.320 


0.320 


0.110 


0.320 


1.30 


0.87 


205 

• 

823 

• 

1109 

• 

65 

• 6 

65 

.6 

0 

.46 

0 

.3 

0 

.02 

0 

.3 

One  of  the  first  questions  a system  designer  may  ask  concerns 
the  relative  advantages  of  beams  spread  either  vertically  or  hor- 
izontally. Since  the  diversity  advantage  stems  from  the  lack  of 
correlation  between  the  two  received  diversity  signals,  the  cor- 
relation fall-off  as  a function  of  squint  angle  in  the  vertical 
and  horizontal  directions  is  the  determining  factor  in  assessing 
the  performance  difference.  The  power  loss  fall-off  as  a ftinction 
of  squint  angle  is  of  the  same  order  for  the  two  directions  and 
therefore  is  not  significant  in  the  selection  process.  The  cor- 
relation fall-off  as  a function  of  squint  angle  can  be 
related  to  the  correlation  distance  in  the  plane  of  the  receive 
antenna.  The  correlation  distance  r^  may  be  defined  in  terms  of 
two  point  antennas,  one  located  on  the  great  circle  path  and  one 
at  a distcuice  r^,  where  the  received  processes  on  the  two  point 
antennas  have  a normalized  correlation  coefficient  of  1/e.  Clearly 
if  the  correlation  distance  is  much  greater  than  the  troposcatter 
system  parabolic  reflector  then  the  correlation  between  squinted 
beeuns  will  be  very  high  whereas  if  the  correlation  distance  is 
small  compared  to  the  parabolic  reflector  the  squinted  beaun  cor- 
relation will  be  small.  The  correlation  distance  is  defined  in 
terms  of  a vector  u in  the  plane  of  the  parabolic  reflector  as 

qp(u)  =*  E{r)  E (r  + u) 

where  E(r,)  is  the  received  electric  field  at  a distance  £ from  the 
center  of  the  scattering  volume.  The  results  of  2.4.1  and  2.5.1 
establish  that  the  above  spatial  correlation  function  is  related 
to  the  scattering  vol\ime  "size”  through  a three  dimensional  Fourier 
integral  u>r' 

cp(u)  « C j'j'J  V(r ' ) e ^ d^r'  . 
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The  integration  is  over  the  scattering  volume  defined  by 

V(r')  which  is  determined  by  the  antenna  gain  patterns  and  the 

refractive  index  spectrum  fall-off  as  a function  of  scattering 

angle.  The  constant  k is  equal  to  2TTf  /c  where  f the  carrier 

o o 

frequency  and  c is  the  speed  of  light.  Now  Fourier  variables  in 
one  domain  are  inversely  related  to  Fourier  variables  in  the 
other  domain,  e.g.,  a short  time  pulse  corresponds  to  a wide 
frequency  band.  Thus  the  correlation  distance  in  a par- 
ticular direction  at  the  receiver  is  inversely  related 
to  the  common  volume  size  in  that  same  direction.  In  Fig. 
3.1  we  ex2Lmine  the  common  volume  dimensions  in  the  vertical  (side 
view)  and  horizontal  (top  view)  directions.  For  angle  diversity 
systems,  the  transmit  beaunwidth  is  typically  smaller  than  the 
minimum  scattering  angle  in  order  to  minimize  the  loss  associat- 
ed with  the  squinted  beam.  The  common  volume  "size"  in  the  hori- 
zontal direction  is  then  limited  by  the  transmit  beeunwidth  and  is 
on  the  order  of  Y^r  = Dr.  In  the  vertical  direction,  the  common 
volume  "size"  is  not  limited  by  the  beamwidth  as  much  as  the  re- 
fractive index  spectrum  fall-off  which  has  a dependence  of  0 
where  0 is  the  scattering  angle.  Since  0 > 0^>  n the  vertical 
common  volume  size  Y^r  is  larger  than  the  horizontal  common  volume 

size  Y r.  The  area  of  useful  scattering  energy  returned  to  the 
H 

receiver  is  figuratively  shown  in  Fig.  3.1  as  a shaded  area.  Since 

Y > Y_  it  follows  that  the  vertical  correlation  distance  will  be 
V H 

less  than  the  horizontal  correlation  distance  for  narrow  beamwidth 
antennas.  Thus  the  vertical  angle  diversity  system  will  result  in 
a lower  correlation  between  squinted  beeuns  than  a horizontal  angle 
diversity  system. 

It  turns  out  that  as  the  transmit  beeunwidth  is  increased 

such  that  n » 0 , the  situation  reverses  and  Y,.  becomes  smaller 
o V 


j 

1 

i 


i 


J 
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than  Y„  with  the  result  of  a smaller  horizontal  correlation  dis- 
tance  than  vertical.  This  result  is  consistent  with  "fat"  trans- 
mitter beeun  measurements  and  analysis  ^3.1]  of  correlation  dis- 
tances in  the  vertical  and  horizontal  planes.  However,  in  any 
practical  application,  the  transmit  beauns  must  be  narrow  in  order 
to  provide  the  required  antenna  gain  for  successful  operation. 

This  qualitative  discussion  of  correlation  distances  is  sup- 
ported by  a quantitative  evaluation  of  the  horizontal  and  vertical 
correlation  distances  in  Section  2.5.4. 

Because  the  vertical  correlation  distance  is  smaller  than  the 
horizontal  correlation  dist?r>^3  for  practical  systems,  the  angle 
diversity  design  should  utilize  a vertical  squinted  beaun  as  the 
first  additional  angle  of  arrival  signal.  Since  a dual  feedhorn 
vertical  angle  diversity  (VAD)  system  increases  the  conventional 
quadruple  diversity  system  to  an  eight  order  system,  the  use  of 
more  angle  of  arrival  signals  is  not  advantageous  both  because  of 
system  complexity  and  the  diminishing  return  from  additional 
diversity.  Thus  we  fix  the  system  design  as  a dual  feedhorn  VAD 
system  and  tvurn  our  attention  to  the  choice  of  vertical  squint 
angle  and  antenna  bores ight/horizon  angle. 

The  angle  between  the  centerlines  of  the  antenna  patterns  of 
the  two  feedhorns  in  a dual  VAD  system  is  defined  as  the  squint 
angle.  When  the  squint  angle  is  appreciably  more  than  a beamwidth, 
the  correlation  between  the  received  angle  diversity  signals  is 
small  but  the  relative  signal  loss  of  the  elevated  beam  may  be 
excessive  due  to  the  increased  scattering  angle.  Decreasing  the 
squint  angle  reduces  this  loss  but  increases  the  correlation  be- 
tween diversity  signals.  The  diversity  combining  loss  due  to 
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signal  correlation  alone  can  be  expressed  for  flat  fading  channels 
at  large  slgnal-to^nolse  ratios  as 

L - -5  LOGd-P^)  (3.3) 

where  P Is  the  normalized  correlation  coefficient  between  two  com- 
plex Gaussian  processes.  This  function  Is  plotted  In  Fig.  3.2. 

The  optimum  squint  angle  for  the  three  example  systems  has  been 
computed  by  using  this  relation  In  conjunction  with  the  calcula- 
tion of  the  correlation  coefflcl^t  and  squint  loss  as  a function 
of  squint  angle.  This  calculation  used  the  turbulent  scatter  model 
and  common  volume  integration  technique  described  In  Section  2. 

The  results  are  shown  In  Figs.  3.3  - 3.5. 

In  confuting  the  diversity  coznblnlng  loss  DCL,  the  squint 
loss  results  In  an  effective  loss  of  S^/2  because  performance 
at  large  slgnal-to-nolse  ratio  SISIR  Is  determined  by  the  geometric 
mean  of  the  diversity  branch  SNR  values  or  equivalently  by  the 
mean  of  the  SNR  values  In  dB. 

In  these  exanples  and  all  others  confuted,  the  optimmn  squint 
angle  was  determined  to  be  approximately  1 beamwldth.  Also  note 
that  the  optimum  Is  relatively  broad  such  that  the  dB  loss  Is  small 
if  the  squint  amgle  should  be  larger  than  one  beamwldth.  Since 
feedhorn  design  Is  more  difficult  with  smaller  squint  angle  re- 
quirement, this  result  has  In^rtant  practical  inpllcations . The 
feedhorn  design  for  the  RADC  test  link  resulted  in  the  measured 
squint  angles  given  in  Table  3.2.  Vertical  polarization  Is  used 
for  signal  reception  In  the  raoc  tests.  The  additional  loss  for 
this  link  resulting  from  a squint  angle  larger  than  1 beamwldth 
can  be  determined  from  Fig.  3.3  to  vary  between  0.1  and  0.4  dB. 
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Fig.  3.4  Optimvim  Squint  Angle,  S.l^epesi 
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Optimum  Squint  Angle,  Oslo-Krlstlansand 


Table  3.2 

Test  Link  Paraoieters 


The  next  question  to  be  addressed  is  the  choice  of  antenna 
pointing  angle,  i.e.,  the  angle  between  the  main  becun  antenna 
pattern  centerline  (boresight)  and  the  horizon.  As  this  angle 
is  decreased,  the  relative  loss  of  the  elevated  beeun  is  reduced 
but  more  of  the  main  beam  pattern  is  blocked  by  the  ground  such 
that  an  optimum  pointing  angle  exists.  Since  operational  angle 
diversity  ^stem  are  utilized  in  both  directions,  the  boresight/ 
horizon  angles  at  transmitter  and  receiver  are  constrained  to  be 
equal  in  the  optimization  search.  Using  the  optimiim  squint  angle 
of  one  beamwidth,  the  main  beam  loss  for  a conventional  system 
was  determined  as  a function  of  antenna  pointing  angle  and  in 
addition  the  diversity  combining  loss  due  to  squint  loss  auid 
correlation  was  computed  from  the  model  of  Section  2.  The  sum 
of  these  losses  establishes  the  optimum  pointing  angles  for  the 
three  example  systems.  The  results  are  shown  in  Figs.  3.6  - 3.8. 

In  these  and  other  examples  the  optimum  boresight/horizon  angle 
for  the  transmit  and  receiver  antennas  was  found  to  fall  between 
1/4  and  1/2  beamwidth  elevation  above  the  local  horizon.  This  result  was 
only  slightly  influenced  by  the  frequency  of  operation  since  lower 
frequency  systems  use  larger  antennas.  It  is  also  worth  noting 
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Pig.  3.7  Antenna  Pointing  Angle  Optimization 


Fig.  3.8  Antenna  Pointing  Angle  Optimization,  Oslo-Kristiansand 


that  the  optimum  pointing  amgle  for  conventional  systems  would 
be  about  1/2  beeunwidth  (main  beam  loss  curves)  but  that  little 
loss  is  incurred  if  the  antennas  aure  not  realigned  for  an  angle 
diversity  application. 


A final  consideration  in  choice  of  antenna  pointing  au:\gle 
is  the  effect  of  increased  multipath  in  both  beeuns  due  to  larger 
pointing  auigles.  Because  digital  modems  can  utilize  the  multi- 
path  to  increase  the  effective  diversity  order  if  the  multipath 
spread  is  not  too  large,  it  is  advantageous  to  choose  a somewhat 
larger  pointing  angle  if  the  nominal  multipath  spread  is  small 
compared  to  the  data  symbol  interval.  However,  for  L-band  systems 
the  multipath  spread,  particularly  in  the  elevated  beam,  may 
exceed  the  multipath  capabilities  of  the  digital  modem.  In  this 
situation  degraded  performance  due  to  intersymbol  interference 
results. 

The  results  from  the  numerical  integration  of  the  common 
volume  using  a turbulent  scattering  hypothesis  are  shown  in 
Figs.  3.9  - 3.11  for  the  example  systems.  The  2a  multipath 
spread  of  the  product  beam  refers  to  the  cross  channel  multipath 
profile  defined  in  [3.2].  The  C-band  test  link  has  small  multi- 
path  spread*  for  data  rates  of  interest  and  because  of  the  small 
dB  loss  associated  with  larger  pointing  angles,  the  optimiun  point- 
ing angle  for  a digital  VAD  system  is  probably  very  close  to  the 
optimum  pointing  angle  for  a conventional  2S/2F  analog  troposcat- 
ter  system,  i.e.,  approximately  1/2  beamwidth.  For  the  L-band 

* Measured  results  from  C-band  systems  tend  to  be  larger  [3.3] 
due  to  the  Inclusion  of  the  terminal  equipment  filter  cheuracter- 
istlcs  amd  measurement  noise. 
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Fig.  3.9  Multipath  Characteristics,  RADC  Test  Link 
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Pig.  3.10  Multipath  Characteristics 
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systems  there  is  considerably  more  multipath  emd  a steeper  slope 
with  pointing  angle.  Optimum  performemce  of  a digital  VAD  sys- 
tem at  these  frequencies  may  result  from  pointing  auigles  some- 
what less  than  1/4  beamwidth. 

In  svimmary  the  optimum  VAD  ststem  utilizes  a squint  emgle 
of  approximately  1 beaunwidth  with  a bores ight/horizon  amgle  ap- 
proximately equal  to  1/2  beeunwidth  for  C-band  digital  systems 
and  somewhat  less  than  1/4  beaunwidth  for  L-band  systems.  In 
all  cases,  however,  the  optimum  configuration  results  in  less 
received  signal  power  in  the  elevated  beeun.  The  use  of  pointing 
angles  with  equal  beampowers. which  has  been  suggested  [3. 4]. is 
suboptimum  because  too  much  power  is  blocked  in  the  main  beam 
by  the  foreground. 

3.2  Long  Term  Variability 

Tropospheric  scatter  systems  are  subject  to  two  fading 
phenomena  ^ short  term  multipath  fading  and  long  term  power 
fading.  The  short  term  fading  of  the  instantaneous  received 
power  within  periods  of  time  ranging  from  less  than  a second  to  many 
minutes  results  from  random  fluctuations  in  the  relative  phasing 
between  component  waves  arriving  at  the  receiver  over  slightly 
different  propagation  paths.  The  long  term  power  fading  results 
from  slow  changes  in  average  atmospheric  refraction,  in  the 
intensity  of  refractive  index  turbulence,  and  in  the  degree 
of  atm^sherlc  stratification.  The  power  fading  is  characterized 
by  hourly  or  diurnal  variations.  The  evaluation  of  troposcatter 
system  performance  is  accomplished  in  part  by  determining  the 
hourly  median  path  loss  where  the  median  is  confuted  to  include 
the  short  term  multipath  fading  and  in  turn  considering  the 
median  path  loss  as  a random  variable  subject  to  a power  fading 
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distribution.  The  median  path  loss  calculation  utilizes  fixed 
values  of  the  mean  and  variance  of  the  refractive  index  and  an 
assunption  on  the  degree  of  atmospheric  stratification.  In 
Section  2 we  have  presented  a prediction  method  for  the  confuta- 
tion of  mediaui  path  loss.  This  method  has  the  following  infor- 
tant  characteristics: 

* The  average  atmospheric  refraction  is  fixed  by 
utilizing  an  effective  earth's  radius  of  4/3  the 
actual  radius  to  account  for  the  mean  refractive 
index. 

* The  intensity  of  refractive  index  tvirbuleurice  is 
fixed  by  the  refractive  index  varizuice  which  is 
chosen  to  correspond  to  dry  winter  afternoons. 

This  period  of  time  generally  experiences  the 
poorest  propagation  conditions. 

• The  atmospheric  structure  is  derived  from  turbulent 
scattering  theory  which  leads  to  a refractive  index 
spectrum  slope  of  n*ll/3.  This  structiire  is  more 
applicable  for  higher  frequency  (>  1 GHz)  tropo- 
scatter  systems  than  the  stratified  layer  assumption 
used  in  the  NBS  prediction  method. 

• The  apertvire-to-medium  coupling  loss  is  included 
as  an  integral  part  of  the  path  loss  calculations. 

For  this  method  the  power  fading  of  the  median  path  loss 

is  determined  by  veuriations  in  the  effective  earth  radius  fac- 

2 

tor  K,  variations  in  the  refractive  index  veuriance  0 , and 

' n 

changes  in  the  atmospheric  structure  leading  to  other  values  of 

the  refractive  index  spectrum  slope  n.  Given  probability 

2 

density  functions  on  the  peurameters  K,  and  n,  the  computa- 
tion of  the  median  path  loss  long  term  distribution  would  be 


* This  time  period  is  referred  to  as  Time  Block  2 (TB2)  in  the 
NBS  prediction  method. 
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straightforward.  Unfortunately  there  is  little  empirical  data 

available  to  derive  such  densities.  Some  experimental  evidence 

and  analysis  [3.5]  indicate  that  the  refractive  index  variance 
2 

is  the  dominating  factor  in  producing  significant  variations 

in  the  median  path  loss.  Development  of  experimental  data  on 
2 

over  long  periods  of  time  for  different  geographical  areas 

would  provide  a basis  for  predicting  long  term  variability  of 

troposcatter  systems.  At  the  present  time  the  only  method  of 

predicting  this  variability  is  to  use  path  loss  data  taken  from 

existing  systems  and  Integrated  into  the  NBS  variability  model 

[3.6].  Much  of  the  empirical  path  loss  data  has  been  taken  from 

systems  with  operating  carrier  frequencies  belcxtr  1 GHz.  The 

performance  prediction  for  new  troposcatter  systems  operating 

in  the  4 to  5 GHz  frequency  region  may  be  subject  to  large 

errors  as  a result.  However,  in  the  absence  of  empirical  data 
2 

on  either  or  median  path  loss  at  these  new  frequencies,  an 
'extrapolation  of  the  NBS  variability  model  is  the  only  realistic 
engineering  choice.  In  the  next  subsection  we  briefly  review 
this  model  and  describe  the  specific  parameters  integrated  into 
the  SIGNATRON  prediction  computer  progreun  for  computation  of 
long  term  variability.  In  the  subsequent  subsection  the  spe- 
cial issues  which  must  be  addressed  in  the  evaluation  of  the 
long  term  variability  of  an  angle  diversity  system  are  con- 
sidered. 

3.2.1  NBS  Long  Term  Variability  Model 

Considerable  experimental  evidence  suggests  that  the  long 
term  distribution  of  the  hourly  median  path  loss  is  normally 
distributed  in  dB.  If  we  denote  L(p)  as  the  hourly  median 
path  loss  in  dB,  which  is  not  exceeded  p%  of  the  time,  the 
normal  distribution  is  defined  as  follows: 
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(3.1) 


100 


= prob  I median  path  loss  < L (p) | 


L(p) 

where  L(50)  is  the  median  path  loss  computed  using  the  method 
of  Section  2 and  is  the  long  term  standard  deviation.  The 
NBS  Long  Term  Variability  Model  r3.6]  uses  empirical  data  to 
determine  L(10)  and  L(90)  from  which  normal  probability  graph 
paper  can  be  used  to  plot  L(p)  at  other  values  of  p.  This 
calculation  includes  the  effect  of  prediction  uncertainty  through 
a parameter  called  service  probability.  The  service  probability 
is  the  probability  that  a new  system  will  meet  the  long  term 
performance  predictions . 

A parameter  called  effective  distance,  d^,  has  been  found 
to  be  superior  to  other  parzuneters  such  as  path  length, angular 
distance, and  distance  between  actual  or  theoretical  horizons, 
in  predicting  the  long  term  variability.  The  effective  dis- 
tance is  defined  in  terms  of  dgj^  which  is  the  distance  for 
equal  diffraction  and  scatter  path  loss  and  dg^  which  is  the 
smooth  earth  distance  between  radio  horizons.  These  parameters 
are  computed  from  the  following  relations 

d - 65(100/f)^^^  hm  (3-2) 

®1 

d - d - 3 -3  -/TIT  (3.3) 

"o  '^e  e 

where  the  effective  antenna  heights  h^  and  h^  are  expressed 
In  meters,  the  path  length  d is  in  kilometer s,an9  the  radio 


[ exp  [(?-L(50))^  / ] d5  . 
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frequency  £ is  in  MHz.  The  effective  disteuice  is  defined  in 
the  NBS  nusdel  as 


! 


130/  [l+(d  - d ) /d]  tan,  d < d 

I «o" 

1 130  + d-d  tan,  d>d 

L «1  ®1 


(3.4) 


A variability  function  y(p)  which  depends  on  the  effective 
distance  and  can  be  corrected  for  frequency  and  climate  effects 
is  used  to  determine  L(10)  and  L(9(J  from  the  relation 


L(p)  - L(50)  - Y(p)  (3.5) 

The  variability  function  is  defined  by 
Y(P)  “ Y^(p,  d^)  g(p,  f) 

where  g(p,  f)  represents  an  average  of  many  effects  that  are 
frequency  sensitive.  The  NBS  varied>ility  model  provides  empiri- 
cal curves  and  an  analytic  function  representation  for  Y^(p,d^). 
In  addition, the  model  uses  a parauneter  V(50,d^)  to  adjust  the 
long  term  reference  median  path  loss  computed  by  the  NBS  method 
to  the  median  loss  L(50)  in  (3.5).  For  dry  winter  afternoons 
this  adjustment  is  0 dB. 

The  analytic  function  representation  for  Y^(p,d^)  is  given 
for  dry  winter  afternoons  (Time  Block  2)  as 

■ [=1  - 'j'V]  “P(-V. 


(3.6) 


Y^dO,  d^) 

O B 


-Yi90,  d^) 

O B 


Y(O.Ol)  - 3.33Y(10)  Y(99.99)  - 2.90Y(90) 

Y(O.l)  ■ 2.73Y(10)  Y(99.9)  = 2.4lY{90) 

Y(l.O)  - 2.00Y(10)  Y(99)  - 1.82Y(90). 

ThM«  calculations  and  the  loss  distribution  resulting  from 
(3.5)  provide  the  long  tern  variability  of  the  median  path  loss 

I 

i 

i 
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for  a continental  temperate  climate  during  winter  afternoons. 
The  NBS  variability  model  provides  numerous  curves  auid  other 
variability  functions  to  accommodate  other  portions  of  the 
season  and  other  climatic  areas. 

The  path  loss  distribution  (3.5)  corresponds  to  a service 
probability  of  50%,  i.e.,  50%  of  the  systems  built  would  exceed 
the  performance  predictions.  Conservative  engineering  practice 
would  recommend  the  selection  of  a service  probability  of  95%. 
For  this  choice  the  loss  distribution  must  be  adjusted  by  the 
prediction  error  according  to  the  NBS  formula 

Lo  95(p)  = L(p)  + 1.65  / 12.73  + 0.12Y^(p).  (3.9) 


The  variability  predictions  used  in  this  study  utilize  the 
NBS  long  term  variability  model  with  a service  probability  of 
95%.  The  median  path  loss  L(50)  is  computed  using  the  method  of 
Section  2, which  corresponds  to  the  loss  for  turbulent  scattering 
conditions  during  dry  winter  afternoons  in  a temperate  continental 
climate.  The  formulas  3.1  through  3.9  are  then  used  to  derive 
the  loss  distribution.  The  mean  distribution  for  evalua- 

tion long  term  variations  in  digital  system  performance  is  de- 
termined from  the  path  loss  distribution  by  the  formula 


^(P) 

o 


G„-  L .^(p)  -^10  LOG  (R.  ) +NF  - 174*^+ 1.6 
T T R .95  i_  D «• 

(3.10) 


where 


a transmit  power  in  dBm 

G.  * transmit  antenna  gain  in  dB 

T 

G * receiver  antenna  gain  in  dB 


i 
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=s  path  loss  not  exceeding  p%  of  the  time  for 
unit  gain  antennas  and  a service  probability 
of  0.95 

» data  rate  in  b/s 

= receiver  noise  figure  in  dB 

a received  noise  power  in  1 Hz  bandwidth 

s factor  relating  median  to  mean  for  a con^>lex 
Gaussian  scatter  channel. 

3.2.2  Long  Term  Variability  in  an  Angle  Diversity  System 

In  space  amd  frequency  diversity  systems,  the  common  scat- 
tering volume  is  virtually  the  same  for  each  diversity  and  hence 
the  long  term  medieui  path  loss  varies  the  same  for  each  diversity. 
In  angle  diversity  system,  however,  the  common  volumes  are  sepa- 
rated and  the  long  term  vsurieOsility  is  not  identical  for  the 
angle  diversity  beams.  The  effect  of  this  decorrelation  of  the 
diversity  power  fading  may  improve  the  system  availability  be- 
cause a power  fade  in  the  main  beam  diversity  is  not  always  ac- 
companied by  a power  fade  in  the  squnited  diversity  beam. Phys- 
ically one  can  imagine  this  situation  in  n vertical  angle 
diversity  system  where  the  inhomogeneous  structure  of  the  at- 
mosphere results  in  say  a larger  refractive  index  variance  or 
superior  atmospheric  stratification  at  heights  corresponding 
to  the  elevated  beam  than  at  heights  defined  by  the  main-beaun 
common  volume.  Since  the  common  volumes  in  a vertical  angle 
diversity  system  are  separated  by  approximately  one  beamwidth 
0 at  a distance  d/2  from  the  link  terminals , this  decorrelation 
results  from  atmospheric  variations  over  distances  on  the  order 
of  nd/2.  For  typical  troposcatter  applications  this  distamce 
is  on  the  order  of  one  or  more  miles  at  a height  above  the 


* Apertxire  to  medixxm  coupling  loss  of  the  actual  diameter 
antennas  is  included  in  this  parameter. 


\95<P' 

NF 

-174dBm 
1.6  dB 
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earth's  surface  also  on  the  order  of  a few  miles.  Experimental 
evidence  verifying  the  long  term  power  fading  decorrelation  in 
angle  diversity  systems  has  been  reported  on  by  Monsen  [3.7]  emd 
Troitskiy  [3.8].  Data  from  these  angle  diversity  experiments  will 
be  used  to  estcdslish  the  availability  improvement  from  power  fad- 
ing decorrelation. 

The  long  term  median  path  loss  can  be  described  by  a multi- 
variated  normal  distribution  when  each  path  loss  is  given  in  dB. 
We  consider  a dual  vertical  angle  diversity  system  with  or  with- 
out additional  frequency  or  space  diversity.  Let  x^  represent 
the  median  path  loss  in  dB  for  diversity  i and  number  of  diversi- 
ties such  that  the  main  beam  diversities  are  odd  and  the  elevated 
beam  diversities  are  even, e.g. , for  a 2S/2F/2A  system  we  have 
x^,  x^  mainbeam  space  diversity  channels 
x^,  x^  mainbeam  frequency  diversity  channels 
Xj,  x^  elevated  beam  space  diversity  channels 
Xg,  Xg  elevated  beam  frequency  diversity  channels. 

Because  of  the  power  fading  correlation  in  space  and  fre- 
quency diversity,  one  has 

X St  X * X 3 X 

1 3 5 7 


and  x^  euid  x^  are  described  by  the  joint  normal  density  function 


f(Xj^,  Xj) 


-x'l^/  2<7 


2m  |r| 


(3.11) 


where 

^ _ r ^1  ” “l 
L *2  " ®2 


“l  at 
aP  a" 

w 
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a mean  value  of  median  path  loss 
for  diversity  1 

m^  * mean  value  of  median  path  loss 
for  diversity  2 

a mainbeam  path  loss  standard 

deviation 

aa  elevated  beam  path  loss  standard 

deviation 

p » power  fading  correlation  coefficient. 

The  NBS  variability  model  discussed  in  the  last  subsection 
provides  a method  for  determining  the  path  loss  standard  devia- 
tion a.  Since  p^o  and  a«l  for  space  and  frequency  diversity 
systems*  determination  of  a completely  specifies  the  long  term 
variability  for  these  systems.  In  angle  diversity  systems  the 
parameters  p and  a will  also  influence  the  long  term  variability 
of  the  effective  received  signal  power. 

In  order  to  assess  ehe  effect  of  decorrelated  power  fading 
and  an  increase  in  the  elevated  beam  standard  deviation,  we  use 
the  mean  diversity  path  loss  in  dB  as  a system  measure  of  per- 
formance. This  measure  corresponds  to  the^  geometric  mean  of 
the  diversity  signal-to-noise  ratios  which  asymptotically  for 
large  signal-to-noise  ratio  governs  the  short  term  performance 
of  both  analog  and  digital  systems.  Thus  we  define  the  mean 
diversity  path  loss 

X - Xj)  (3.12) 

and  a system  availability  fttnction  as  x(p),  where  x(p)  is  the 
mean  diversity  path  loss  not  exceeded  p}(,  of  the  time,  viz., 

I 
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100 


prob  (mean  diversity  path  loss  < x(p)| 


Since  the  stun  of  two  normal  random  variables  is  also  normal* 
we  have 


x(p) 

100 V"'!'  ^ 


(3.13) 


where  the  long  term  mean  and  variance  are  equal  to 

X = I (*1  + *2^ 

“x^  ■ [ '*1  - ”l>  <*2  - V 

2 

= (1  + 2Pa  + a ) . 

4 

The  percent  availability  for  the  angle  diversity  system  is 


100 


1 -Q 


x(p)  - (mj^+m2)/2'' 


Lra 


l+2Pa+a 


where  Q(*)  is  the  normal  distribution  function 

1 P 
0(u)  « I 


-V 


•J  2Tr 


dv 


u 


(3.14) 


(3.15) 


Equation  (3.14)  for  angle  diversity  should  be  contrasted 
with  the  result  for  space  or  frequency  diversity  when  P «o,  a«l> 
viz. , 
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(3.16) 


[ 


p » 100 


Note  that  for  the  angle  diversity  system  is  not  neces- 
sarily equal  to  for  the  space  or  frequency  diversity  systems. 

As  an  example,  when  angle  diversity  is  used  in  place  of  frequency 
diversity,  the  additional  power  aiii>lifier  required  for  frequency 
diversity  results  in  a 3 dB  improvement  for  angle  diversity,  i.e«, 
m^^  ■ m^  + 3. 

In  (3.14)  one  can  see  that  the  effect  of  the  power  fading 
correlation  coefficient  is  to  improve  system  availability  with 
decreasing  values  of  the  coefficient  but,  on  the  other  hand,  an 
increase  (a  > 1)  in  the  elevated  beam  path  loss  standard  devia- 
tion decreases  the  system  availability.  In  order  to  evaluate 
the  result  of  these  competing  effects,  we  examine  some  previous 
angle  diversity  experimental  data. 

Long  term  variability  data  was  taken  on  the  Bell  Laboratories 
experimental  angle  diversity  link  reported  on  by  Monsen  [3.7]. 

The  important  characteristics  of  this  link  were 

• triple  vertical  angle  diversity 

• 2.17  (Sz  operating  frequency 

• 179  statute  mile  path 

• 28'  transmit  antenna,  50'  steerable  receive  antenna 

• 1 Kw  transmit  power 

• 1.6  bezunwidth  separation  between  beams. 

Measurements  of  hourly  median  signals  once  or  twice  a week 
over  a period  from  September  1970  through  February  1971  provided 
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a set  of  43  independent  hourly  median  seunples.  The  long  term 
★ 

received  power  means,  standard  deviations,  and  correlation  co- 
efficient for  the  lowest  two  vertical  beauns  have  beai  computed 
and  are  given  in  the  following  table  as  a function  of  amtenna 
pointing  angle  relative  to  the  horizon. 

Now  let  us  compare  the  availad>ility  of  a space/ frequency 
diversity  system  with  space/angle  diversity  system.  We  select 
an  antenna  pointing  angle  of  0.5  beaonwidths  as  both  the  results 
of  Section  3.1  amd  the  data  in  Table  3.3  suggest  that  performance 
will  be  maximized.  For  a space/ frequency  diversity  system 

2S/2F  diversity:  m^  = m^  = -86.7  dBm 

a = 1 P = o 


P 


100^1  - q(- 


+ 86.7N"1 
6.1 


For  the  2S/2A  diversity  system,  the  second  power  amplifier 
used  to  derive  the  other  frequency  in  the  2S/2F  system  can  be 
tuned  to  the  same  frequency  with  a resulting  increase  in  received 
power  of  3 dB.  Thus  for  the  space/angle  diversity  system 

2S/2A  Diversity  » m^^  * -83.7  dBm 

m^  » -90.0  dBm 

9 » 6.1 

a » 7. 2/6.1  » 1.18 
P - 0.743 

and  from  Eq.  (3.14) 


P 


100  1 


Q 


X + 86.85  \1 
.1.0186x6.1 /j  * 


Received  power  atnd  path  loss  are  used  interchangeably  here  since 
they  are  related  by  an  additive  constant. 


J 
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The  angle  diversity  system  has  slightly  decreased  availability 
because  its  mean  diversity  received  power  is  0.15  dB  less  and  be- 
cause the  effective  standard  deviation  is  increased  by  a factor  of 
1.018.  Thus  in  some  cases  the  squint  loss  and  increased  standard 
deviation  of  the  elevated  beam  path  loss  overcome  the  availability 
improvement  due  to  the  3 dB  power  cunplifier  gain  and  the  power  fad- 
ing decorrelation.  For  this  experimental  system,  however,  the 
squint  loss  was  considerably  larger  than  would  be  designed  for  an 
operational  system. 

When  the  squint  loss  is  equal  to  6 dB  there  is  exact  cancel- 
lation of  the  3 dB  power  amplifier  gain  with  the  mean  diversity 
loss  associated  with  squinted  dual  angle  diverssty  system.  For 
operational  systems  the  net  gain  due  to  these  two  factors  will 
normally  be  positive  since  the  antenna  beams  can  be  placed  closer 
together  than  the  1.6  beamwidth  separation  in  the  BTL  experimental 
system.  The  net  gain  due  to  the  power  fading  correlation  coef- 
ficient and  elevated  beam  standard  deviation  is  more  difficult  to 
assess  because  of  the  rather  limited  empirical  data.  The  BTL 
data  suggests  that  the  two  effects  tend  to  cancel  each  other  out. 
Troitskiy  [3.8]  has  presented  data  on  three  links  which  show  a 
smaller  power  fading  correlation  coefficient  than  found  from  the 
BTL  data.  Table  3.4  stnnmarizes  Troitskiy’s  experimental  results. 
The  power  fading  correlation  coefficient  was  measured  between 
envelope  powers  in  watts  rather  than  dBm  and  is  designated  p. 

Ir  the  Appendix  we  show  that  p is  ovezbounded  by  p but  the  bound 
is  tight  and  to  a good  approximation  p » p. 

It  is  STirprising  that  the  results  in  Table  3.4  show  a smaller 
power  fading  correlation  coefficient  when  these  links  have  the 
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beams  closer  together  in  the  sky  than  the  beams  were  in  the  BTL 
experimental  link.  Also  it  is  vmfortunate  that  Troitskiy  did 
not  report  on  data  from  which  the  elevated  beam  standard  devia- 
tion can  be  .culated.  A major  goal  of  the  AAC  field  test 
program  will  be  to  obtain  estimates  of  both  the  power  fading 
correlation  coefficient  and  the  elevated  beam  standard  deviation. 

This  data  will  be  used  to  help  resolve  what,  if  any,  avail- 
ability improvement  can  be  attributed  to  long  term  decorrelation 
of  angle  diversity  beams. 


Table  3.4 

Experimental  Results,  Ref.  Troitskiy  [3.8] 
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SECTIOtT  4 

DIGITIOi  SYSTEM  PERFORMANCE 

The  bit  error  rate  performance  of  a digital  data  system 
operating  over  an  angle  diversity  troposcatter  link  is  dependent 
on  the  received  power  at  each  feedhorn,  the  correlation  between 
diversity  chamnels,  and  the  multipath  chzuracteristics  of  the 
composite  diversity  system.  As  one  of  the  end  products  of  this 
study,  we  will  develop  the  mathematical  formalism  to  predict  (1) 
the  average  bit  error  rate  as  a function  of  the  hourly  median 
received  power  for  the  main  beam  diversity  channel  and  (2)  the 
long  term  probability  that  the  average  bit  error  rate  is  exceeded. 

The  formalism  will  be  converted  into  computer  progreuns  for  con- 
venient calculation  of  these  results  for  specific  angle  diversity 
proposed  applications . 

The  first  goal  is  concerned  with  evaluating  the  short  term 
(less  than  a minute)  fading  chauracteristics  of  the  troposcatter 
channel.  The  second  goal  deals  with  the  availability  of  the 
digital  troposcatter  system  as  a result  of  hourly  and  dally  vetria- 
tions  in  path  parameters.  The  path  loss  auid  multipath  prediction 
model  developed  under  this  program  provides  the  basis  for  the 
calculation  of  the  short  term  path  parameters  while  the  results 
from  3.2  euid  the  NBS  long  term  variability  data  will  be  used  to 
assess  the  availability  question. 

Confutation  of  bit  error  rate  statistics  Implies  a choice 
of  digital  modulation  format  and  receiver  processing  technique. 

We  have  selected  QPSK  for  the  modulation  format  and  a generalized 
adaptive  decision- feedback  equalizer  (DFE)  for  the  processor. 

QPSK  is  the  present  choice  of  the  two  developed  troposcatter  modems, 

MD-918  and  the  DAR-4,  because  of  its  bandwidth  efficiency,  nearly 


f 


I 


r 


constant  envelope  characteristic,  and  optimal  detection  capability. 
The  OFE  has  been  shown  in  extensive  field  testing  of  the  ND-918 
modem  to  provide  excellent  performance  over  an  extremely  wide 
range  of  channel  conditions  such  as  might  be  anticipated  in  an 
angle  diversity  application.  The  MD-918  DFE  modem  is  also  the 
modem  to  be  used  in  the  AAC  field  tests. 

The  general  problem  of  computing  the  short  term  bit  error 
rate  (BER)  distriubution  for  a fading  multipath  channel  is  a dif- 
ficult if  not  impossible  task.  One  can  determine  the  BER  dis- 
stribution  for  a flat  fading  channel  with  diversity  order  D in  a 
straightforward  manner  [4.1]  but  the  introduction  of  in^licit 
diversity  due  to  multipath  precludes  a closed  form  transforma- 
tion of  variables.  The  short  term  mean  BER  for  a space  or  fre- 
quency diversity  system  can  be  determined  following  the  procedure 
developed  in  [4.2]  with  a resulting  expression  in  the  form  of  a 
determinant.  The  BER  distribution  can  be  estimated  by  finding 
the  mean  BER  and  then  defining  an  equivalent,  non  integer,  flat 
fading  diversity  system  by  matching  the  mean  BER's  in  the  region 
of  interest.  The  effective  order  of  diversity  for  the  equiv*- 
alent  flat  fading  system  could  then  be  used  to  derive  the  short 
term  BER  statistics.  The  effective  diversity  order  would  reflect 
the  correlation  loss  between  diversity  channels,  the  squint  loss 
of  elevated  beams,  implicit  diversity  due  to  multipath,  euid  the 
intersymbol  interference  penalty  which  limits  the  usable  in^licit 
diversity.  Thus,  the  computation  of  mean  bit  error  rate  is  a 
logical  first  step  in  the  evaluation  of  system  performance.  In 
the  subsection  to  follow,  the  analysis  developed  in  [4.2]  is 
extended  to  include  correlated  explicit  diversity  channel  systems 
such  as  angle  diversity. 
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4.1  Mean  Bit  Error  Rate  Derivation 

4.1.1  Communication  System  Definition 

The  communication  system  model  under  consideration  is  shown 
in  Fig.  4.1.  Complex  notation  is  used  to  represent  in-phase  auid 
quadrature  components  and  e3q>licit  modulation/deroodulation  opera- 
tions are  not  shown.  The  data  to  be  trzmsmitted  (s^)  are  selected 
from  the  set  (±l+j)  for  a quaternary  phase- shift-keying  system. 

The  trauismitted  waveform  is 

CD 

a^f  (t-kT)  (4.1) 

where  T is  the  interval  between  successive  symbol  trauismission, 
i.e.,  symbol  interval.  The  data  rate  for  this  QPSK  system  is 
2/T  bits/second.  The  transmitter  impulse  response  f (t)  is  de- 
fined as  a unit  energy  real  function. 

The  fading  multipath  channel  can  be  represented  by  an  ensemble 
of  zero-mean  random  functions  with  complex  Gaussian  statistics. 
Individual  diversity  channels  in  an  angle  diversity  system  will 
exhibit  correlation  between  each  other  which  can  be  completely 
described  by  the  second  moment  of  the  channel  impulse  response  h^  (t) , 

h^(t)  hj(T)  = 6(t-T)  i»  j * 1,2,..,D, 

(4.2) 

where  6 (t)  is  the  Impulse  fhnctlon,  is  the  mean  re- 

ceived energy  per  bit  for  the  main  beam  diversity  channel, 
represents  a relative  received  power  difference  between  the  diversity 
channels,  and  p^^  (t)  is  a unit  area  fttnction  equal  to  the  mean 
squared  envelope  of  the  incise  response.  It  represents  the  aver- 
age cross  power  response  at  each  delay  value  to  impulse  function 
excitation  and  will  be  referred  to  as  the  cross  multipath  profile. 
Twice  the  RMS  width  of  the  profile  is  the  cross  multipath  spread, 
20^^.  For  convenience  the  main  beam  diversity  channel  is  designated 
as  channel  1 which  requires  a typical  angle  diversity 
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Pig*  4.1  Communication  System  Model 


application  utilizing  dual  vertical  feedhorns  on  each  antenna 
with  odd  nvuiibers  designating  main  beam  diversities  and  even 
nunibers  designating  elevated  beam  diversities, the  characteriza- 
tion has  the  special  form 


A.  . = 1 

# 

Pii(t)  * P^(t)  ; 

i 

odd 

(4.3.1) 

A.i  ^ 1 

9 

P^i(t)  = P2(t)  ; 

i 

even 

(4.3.2) 

\,i+l  s 

1 ' 

Pi,i+1  “ ^12^^^  * 

i 

odd 

(4.3.3) 

p. . (t)  = 

0 , 

i j and  i (odd) 

j-1 

(4.3.4) 

Conditions  (4.3.1)  (4.3.4)  state  that  all  main  beam  diversities 
have  equal  statistics  euid  are  independent  of  each  other.  Conditions 
(4.3.2)  2Uid  (4.3.4)  state  that  all  elevated  beauns  have  equal  statis- 
tics and  are  independent  of  each  other.  Finally,  conditions  (4.3.3) 
and  (4.3.4)  state  that  the  cross  multipath  profile  is  the  saune  for 
each  feedhorn  pair  but  that  the  correlation  between  feedhorn  ports 
on  different  antennas  is  zero.  Note  that  the  elevated  beam  in  a 
vertical  angle  diversity  system  generally  has  less  received  power 
because  of  its  laurger  scattering  angle.  The  results  from  Section 
3.1  contend  that  in  a duplex  angle  diversity  system  it  is  not 
advantageous  to  point  the  antennas  down  to  equalize  the  received 
powers  in  the  two  feedhorns. 

The  definition  of  Ej^  as  the  received  energy  per  bit  for  the 
main  beam  diversity,  » 1,  follows  from  the  fact  that  one  re- 
ceived bit  is  represented  by  the  pulse 


+ 1 J f(t-u)  hj^(u)  du 

which  has  energy 

« • • 00 

J dt|  J f(t-u)h  (u)du|^  ■ A,,  J dt  J f^(t-u)  E^^p.  (u)  du 

m 

■ ®b  I Pi'"’ 
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The  enseiDble  representation  is  used  here  because  the  analysis 
to  follow  assumes  a slowly  varying  multipath  channel  which  the 
equalizer  can  track.  Hence  the  mean  BER  performance  Is  determined 
by  computing  the  channel  ensemble  average  of  the  bit  error  pro- 
bability for  the  Decision-Feedback  Equalizer  (DFE) . An  overbar 
Is  used  to  denote  channel  ensemble  averages.  Brackets  <)  will  be 
used  to  denote  average  over  noise  or  source  statistics  for  a par- 
ticular channel  realization.  For  the  white  Gaussian  noise  chan- 
nels the  additive  noise  terms  n. (t)  are  zero  mean  and  have  second 
moments 

<n^(t)  n^*(T)>  - fi(t-T).  1,  j«l,2,...D, 

where  Is  the  noise  spectral  density  In  watts/Hz. 

Each  forward  filter  In  a realizable  DFE  consists  of  a finite 

length  tapped  delay  line  filter  with  tap  spacing  of  r seconds. 

s 

If  there  are  "late"  taps  and  "early  taps",  the  tapped  delay 
line  filter  Impulse  response  Is 

^2 

w. (t)  = r W*  6 (t-kT  ) (4.4) 

^ k-K^  ® 

where  the  tap  gain  value  Is  chosen  with  a complex  conjugate  for 
notational  convenience  later.  The  receiver  filter  has  impulse 
response  f(-t)  which  matches  the  treuismit  pulse  waveform.  For 
this  choice, In  the  absence  of  multipath, and  after  time  synchroni- 
zation, the  optimum  tap  weights  for  the  main  beam  diversity  forward 
filter  reduce  to  one  tap  on  and  the  rest  off.  The  output  of  the 
receiver  filter  is  sampled  at  the  symbol  rate  1/T  to  produce  the 
sample  X . The  backward  filter  correction  sample  la  summed  with 

IQ 

the  sampled  receiver  filter  output  X^  to  produce  the  detection 
variable  z at  the  mth  symbol  epoch.  The  backweurd  filter  weights 

IQ 

previous  decision  by  the  complex  weight  b^  to  form  the  sample 

output 


(4.5) 


'm 


B 

r 

i-l 


b.  s - 
1 m-1 


at  the  decision  time  for  the  s symbol. 

m 


The  decision  on  s is 

m 


denoted  as  s . The  decision  process  teikes  the  form 


m 


z = X + y 
m m '‘m 

S » sgn  (Re(z_,))  + jsgn  (Im(z^)) 
mm  m 

where  sgn  is  the  signvim  function. 

The  parameters  of  the  DFE  are  the  number  of  forward  filter 

taps  K,  the  forward  filter  tap  spacing,  t and  the  number  of  back- 

s 

ward  filter  taps,  B.  The  optimum  DFE  requires  K » * » B and 

T « 1/W  where  W is  the  two-sided  bandwidth.  A practical  choice 
s 

of  parameters  for  troposcatter  channels  was  determined  from  com- 
puter simulation  of  a fading  channel  equalizer  application  to  be 

K » 3 « B and  T » T/2. 

s 

4.1.2  DFE  Performance  Analysis 

Previous  methods  of  determining  average  bit  error  rate  (ABER) 
performance  of  practical  equalizer  struct\ires  have  been  restricted 
to  Monte  Carlo  simulations  [4.3,  4.4]  using  an  ensenible  of  multi- 
path  channels.  An  analytic  approach  for  calculating  the  ABER  is 
con^licated  by  the  presence  of  an  intersymbol  interference  term  in 
the  signal-to-noise  ratio  (SNR)  expression  for  a particular  chan- 
nel realization.  Because  of  this  term,  an  average  over  the  channel 
ensemble  is  a formidable  task.  On  the  other  hand,  omission  of  the 
intersymbol  interference  term  leads  to  a mean  BER  which  for  a fixed 
SNR  monotonically  decreases  with  increasing  multipath  spread.  For 
any  practical  equalizer  for  a fixed  SNR  the  mean  BER  will  initial- 
ly decrease  for  increasing  multipath  spread  but  then  increase  as 
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the  multipath  exceeds  the  equalizer's  capability  to  mitigate  the 
intersymbol  interference.  Elimination  of  the  intersymbol  term 
provides  a convenient  lower  bound  which  both  shows  the  intersymbol 
interference  penalty  and  is  am.  accurate  performance  estimate  for 
small  multipath  spreads.  However,  the  lower  bound  is  too  loose 
for  performance  calculation  when  the  rms  multipath  spread  is  on  the 
order  of  the  forward  filter  width.  The  analytic  procedure  to  be 
developed  allows  accurate  calculation  of  the  lower  bound  for  no 
intersymbol  interference  and  by  an  approximation  of  the  intersym- 
bol interference  effect  provides  a good  performance  estimate  for 
the  large  multipath  case. 

For  notational  simplicity  the  analysis  consideres  a nondiversity 
channel  and  extends  the  final  results  for  higher  order  diversity  sys- 
tems including  angle  diversity.  The  received  signal  has  the  form 

r(t)  » I Sji^  Jh(u)  f(t-mT-u)  du  + n(t)  (4.6) 

m ■— **  — cB 

We  define  the  combined  transmitter  and  receiver  response  as 

OB 

g(t)  « f(t+u)  f(u)  du 
and  its  convolution  with  the  channel  as 


g(t-u)  h{u)  du. 


(4.7) 


For  a sampling  time  t^,  the  detection  variable  z has  the  form 


2 • f — 

!k  It.-!-  ■ 


where  the  X-Kj^+K^+l  forward  filter  taps  are  apportioned  between 
"past"  taps  and  Ki  "future"  taps  and  one  center  tap  . The  noise  variate 
v(t^)  • jn(t^-u)  f(u)  du 


is  zero-mean  Gaussian  with  second  moment 


There  exists  an  obvious  representation  in  a K^Kj^+Kj+l  dimensional 
colvunn  vector  form,  viz.. 


where  the  accent  meurk  refers  to  complex  conjugate  transpose.  The 
vector  w represents  the  forward  filter  complex  tap  gains,  3 rep- 
resents the  sampled  continuous  filter  response  for  the  mth  trans- 
mitted symbol,  and  v is  a zero— mean  Gaussian  noise  vector  process 

with  positive  definite  covariance  matrix  N G with  element  values 

o o 

N^g(kT-iT). 

For  QPSK,  every  other  bit  decision  is  made  on  the  real  part 
of  z and  there  is  quadrature  symmetry  and  thus  the  real  noise  power 
affecting  that  decision  is 

6^  =«  <Re(  w'v)  > = “T"  w*G  w (4.9) 

O 2 •“  or* 

If  one  assumes  a =0,  m ^ 0,  there  is  no  intersymbol  inter- 
m 

ference  and  the  analysis  for  the  mean  BER  leads  to  an  SNR  expression 
for  a particular  chzmnel  realization  which  is  a quadratic  form  of 
the  type  g”^  where  g^  is  a complex  Gaussian  vector  with 
statistics  determined  by  the  channel  and  G^  is  a positive  definite 
matrix.  After  a diagonalization  procedure  the  average  of  the  bit 
error  probability  fiinction  of  the  SNR  can  be  computed  for  bbth 
coherent  and  differentially  coherent  detection.  This  calculation 
leads  to  the  lower  bound  expression.  When  the  inter symbol  inter- 
ference is  present  an  effective  signal-to-noise  ratio  can  be  de- 
fined which  leads  to  a quadratic  form  but  both  the  vector  and  the 
matrix  have  random  components.  A method  of  averaging  over  a func- 
tion of  this  quadratic  form  is  not  apparent.  One  course  worth 
considering  is  to  approximate  the  intersymbol  Interference  effect 
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in  a manner  which  will  modify  the  G matrix  but  keep  its  determin- 

o 

istic  nature.  This  can  be  accomplished  if  the  following  two  ap- 
proximations are  made.. 

(1)  Assume  the  interference  symbols  s , m ^ Q,  axe  Gaussian 
distributed  rather  than  complex  binary. 

(2)  Approximate  the  matrix  contribution  to  G due  to  inter- 
symbol interference  by  its  mean. 

This  approach  assumes  that  the  intersymbol  interference  after 
equalization  in  a fading  multipath  channel  can  be  well  approximated 
by  an  equivalent  non-fading  additive  Gaussian  noise  term.  Since 
most  of  the  bit  errors  are  a result  of  fades  in  the  desired  sig- 
nal, gross  approximation  in  the  interference  characteristics  do 
not  significantly  alter  the  mean  BfSR. 

One  can  reasonably  asstzme  that  the  backward  filter  cancels 
the  interference  contribution  due  to  past  symbols,  i.e.,  symbols 
which  arrive  before  the  syndsol  cvirrently  being  decided  on.  If  we 
arbitrarily  take  the  symbol  selected  for  decision  as  s^,  the  re- 
maining interference  is  classified  as  future  and  is  given  by 

M 

2 ® • (4.10) 

msl 

The  sum  can  be  truncated  after  a few  terms  as  the  dot  product 

will  disappear  for  large  i due  to  the  finite  duration  of  the  pulse 

function  q(t).  The  first  approximation  is  to  take  s.  as  a zero 

^ 2 

mean  complex  Gaussleui  with  mean  magnitude  squared  value  Y . This 
quantity  should  be  less  than  two  as  approximating  a unit  magnitude 
binary  variable  by  a unit  variance  Gaussian  variable  will  certain- 
ly lead  to  pessimistic  performance  results  due  to  the  tails  of  the 

2 

Gaussian  distribution.  A choice  of  Y • 1/2  was  found  to  provide 
excellent  agreement  between  calculated  and  measured  values. 

Equation  (4.8)  can  be  rewritten  as 

» - «■  (sg  ^ * %) 
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and  since  the  additive  noise  contribution  from  v is  independent 
of  the  source  digit  contribution  from  V]^»  the  effective  noise 
power  for  the  decision  on  the  real  part  of  z is 

c2 


6^  ■ <^Re^w'  (v  + v^)^ 


6^  * -2  Q 


N 


W*  G W + ” w'  (v-V-  ' ) w “ ““  w'  Gw 
— cy-  2 — — 1—1  — 2 — — 


(4.12) 


Performing  the  indicated  averaging  with  respect  to  the  source 
digits,  one  finds  from  the  definition  of  v^  in  (4.10)  the  covariance 

M 

(4.13) 


matrix 


<ZlVi'  > . E Si  Si 

m“l 


The  signal-to-noise  ratio  (SNR)  affecting  the  decision  on  the  real 
part  of  z is  then 


26‘ 


(4.14) 


which  we  seek  to  maximize  as  a function  of  the  forward  filter 
weight  vector,  w.  consider  the  generalized  dot  product  defined 
on  the  positive  definite  matrix  G,  l.e.. 

(«»  Z)  ® a’Gv.  (4.15) 


Equation  (4.14)  can  be  written  as 


P ■ 


CRe(w,G"^g^) 
(w#i£)* 


(4.16) 


and  by  a generalization  of  the  Schweurtz  Inequality  we  obtain 

< L ' ^ L te-i 

“ N-  (w,w)  - N 


with  equality  if  and  only  if 


So’ 


(4.17) 
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-1 

w ^ * G q 
—opt  -Sc 


(4.18) 


Note  (4.18)  results  in  a pure  real  term  in  (4.16)  to  satisfy 
the  Re  ( • ) requirement . 


The  maximum  SNR  is  then 


fb  • -1 

^m2uc  * N ^ ® ^ 

o 


(4.19) 


As  a check  on  the  normalization,  if  a ■ 0 and  t «■  0 we  have 

o 


from  Eq.  (4.7) 


q(-kT)  a q(-kT)  * G 


(4.20) 


and  then 


0 = ^ g(0)  . ^ 

max  N ' N 
o o 


(4.21) 


as  is  expected  in  the  absence  of  multipath. 

We  will  now  establish  that  the  mean  square  error  function 


<e^)  » ^ < Iz  - ds 1^> 


(4.22) 


which  is  minimized  by  the  decision- feedback  equalizer  leads  to  the 
same  set  of  forward  filter  weights  as  maocimization  of  the  SNR  in 
(4.16)  for  the  appropriate  choice  of  scale  factor  d.  After  squar- 
ing amd  averaging, equation  (4.22)  becomes 

(e^>  ■ Ej^w'q^q^'w  - 2 d Re^w'q^^+'^^  w'  + d^ 


2 

w*  q^-d  1 + N^  w'  Gw. 


(4.23) 


Minimization  of  (4.23)  is  equivalent  to  the  LaGrange  multiplier 
problem 

min  I” w'Gw  - 28  w’  gjj  (4.24) 


which  has  solution 


w - 8 G 
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(4.25) 


Since  Eq.(4.18)  and  (4.25)  are  the  same  except  for  etn  unimportant 
scale  factor,  we  have  established  that  the  maximum  signal-to-noise 
ratio  for  a minimvim  mean-square  error  DFE  and  a particular  channel 
realization  is  the  quadratic  form  given  by  Eq.(4.19).  The  matrix 
G is  defined  from  (4.12)  and  (4.13)  as 


\ 2 
O msal 


M 


(4.26) 


where  G is  a function  of  the  trcuasmitter  impulse  response  and  is 
o 

a deterministic  matrix.  The  second  term  in  (26)  is  a random  matrix 
with  statistics  determined  by  the  channel  ensemble.  Using  G in 
this  form  for  calculation  of  the  mean  BER  requires  an  average  over  a 
function  of  a quadratic  form  which  has  both  random  vectors  and 
random  matrix.  The  second  approximation  used  to  compute  this 
average  is  to  replace  the  random  interference  term  by  its  mean, 

E.  , M T 

<3  i a + ^ t SL%  (4.27) 

o ma*! 


The  matrix  G is  now  a deterministic  function  of  the  transmitter  emd 
receiver  filter  characteristics  and  the  channel  statistics.  The 
vector  is  a random  vector  associated  with  the  fading  channel 
ensemble.  The  average  bit  error  rate  can  now  be  computed  by  taking 
the  average  over  the  channel  ensemble.  The  bit  error  probability 
for  QPSK  for  any  mendber  of  the  channel  ensemble  is  from  (4.11)  and  (4.19) 


h pr  C?  > 1}  + *5  pr  {?  < -1} 


(4.28) 


where  ^ is  a zero-mean  Gaussian  random  variable  with  standard  devi- 


ation equal  to  Jl/lP 


The  bit  error  probability  is 


^ exp  (-u^/2)  du  ■ erfc 


(4.29) 
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The  STibscript  c denotes  that  the  detection  process  was  co- 
herent, i.e.  the  receiver  knows  the  transmitted  carrier  phase. 

Differential  detection  of  phase-shift-keyed  signals  yields  a bit 
error  probability  [4.5] 

-P 

Pd  - ^5  e (4.30) 

The  tap  gains  of  an  adaptive  MMSE  equalizer  remove  the  phase 
and  frequency  difference  between  the  transmitter  and  receiver 
carrier  clocks  to  an  ambiguity  of  180°  for  binary  transmission 
and  90°  for  quaternary  transmission.  This  ambiguity  is  eliminated 
by  differentially  encoding  the  transmitted  data.  The  performance 
of  coherently  detected  and  differentially  encoded  data  signals  for 
a fixed  is  approximately  double  the  error  rate  given  by (4. 29). 

When  practical  degradations  are  included  in  the  analysis  the  per- 
formance is  closer  to  the  differential  detection  expression  (4.30). 

Thus  the  ensemble  average  of  p^  is  of  Interest  as  a fading  channel 
performance  bound  and  the  ensemble  average  of  p^  is  of  interest  as 
an  estimate  of  realizable  modem  performance. 

Let  y(x)  be  the  probability  density  function  for  p and 

nisx 

Y(S)  its  Laplace  transform.  We  want  to  calculate  the  ensemble 
averages 

CB 

P_  ■ F *5  erfc  -Jx  y(x)dx  (4.31) 

C vO 

CB 

^d  ■ Jo  **  *”*  y(x)<ax  - % Y(i)  (4.32) 

The  Laplace  transform  Y(s)  is  easily  obtained  after  a diagon- 
alization  of  the  quadratic  form  in  Eq.  (4.19). 

* ’ ^ So  - ^ 

Where  a is  a zero-mean  complex  Gaussian  vector  with  a diagonal 
covariance  matrix  Tij  ■ The  diagonalization  resulted  from  | 

the  transformation  1 


J 
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> 


a = MG  ^ q 
— -ao 

and  M is  the  normalized  orthogonal  matrix  for  the  matrix 

g"**  i.e. 

g“^  q q'  g”^  M = 


M'M  « 1. 


(4.34) 

(4.35) 

(4.36) 

(4.37) 


Thus  the  eigenvalues  axe  also  the  eigenvalues  of  the  sym- 
metric matrix  G“^  q q ‘ G“^  and  the  xinsymmetric  matrix  G“^  q q' . 


Since  the  con^nents  of  a are  uncorrelated  Gaussian,  the  probabil- 
ity density  for  is  the  convolution  of  exponential  densities 

which  yields  a Laplace  transform  in  product  form. 

Kj  E. 

y(s)  - fT  (i  . - X ,)• 

O 

Tc— K, 


(4.38) 


The  meam  BER  for  OPSK  is 


*'2  E -1 

5^  - i5Y(i)  • Ij  n (i  + 5C  \) 


(4.39) 


k«-K, 


For  coherent  detection  it  is  necessary  to  find  the  coeffi- 
cients A.,  in  a partial  fraction  e3q>emsion  of  Y(s),  i.e., 
i 

^2  E.  -1 

Y(8)  - t A.  (l  + — X s)  . 

Ic— ^ No  k / 

A rectirsive  method  for  finding  the  coefficients  and  the  resulting 
error  probability  for  coherent  detection  is  detailed  in  [4.4]  where 
the  mean  BER  for  the  infinite  length  optimum  single  pulse  receiver 
was  determined.  We  will  concentrate  on  the  DPSK  expression  (4.39) 
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because  the  DPSK  bit  error  probability  curve  for  the  Gaussian 
noise  non  fading  channel  closely  matches  the  HER  for  the  ND-918 
modem,  i.e.,  the  implementation  loss  and  differential  encoding 
loss  in  the  CPSK  detection  system  is  approximately  equal  to  the 
loss  incurred  in  ideal  DPSK  relative  to  ideal  CPSK.  Figure  4.2 
illustrates  this  result.  Since  in  a fading  channel  application 
most  of  the  bit  errors  will  occur  during  periods  of  low  SNR,  the 
correspondence  between  ideal  DPSK  and  the  measured  DFE  modem 
should  be  closest  at  low  SNR. 


For  a DPSK  system  or  for  the  use  of  the  DPSK  es^ression  to 
approximate  the  Gaussian  noise  non  fading  performzuice  of  a practi- 
cal modem,  calculation  of  the  average  bit  error  probad^ility  does 
not  require  the  determination  of  the  eigenvalues.  Note  that  (4.39) 
has  the  determinant  form 


p^  » h(dt|l  + ^ g'^cI)"^ 


(4.40) 


where  the  covariance  matrix  C is  a function  of  the  transmitter 
filter  multipath  profile,  and  the  sampling  time  t . 


q(t  -kT 
o s 


* 


-^•T  ) 

O S 


-J 


g(t  -kT  -u)g(t  —tT  -u)p(u)du. 

OS  OS 


(4.41) 


The  matrix  G as  defined  in  (4.27)  has  the  form 


V-  * N ”,  St <%-*'>•  <*•«> 

o m»l 


The  matrices  (4.47)  and  (4.48)  and  the  signal- to-noise  ratio  com- 
pletely determine  the  mean  BER.  The  extension  to  a diversity  con- 
figuration with  a forward  filter  associated  with  each  diversity 
input  and  zurbitrary  correlation  between  diversity  channels  Involves 
a straightforward  generalization  of  the  matrices  C and  O.  Let 


BIT  ERROR 


c ® 


■c<“>  c<“'  ■ 

,,(21)  ,,(22)  __  ,,2D) 

g(21),,(22)  __g(2D) 

e 

A 

G = 

• 

• 

f 

• 

^(Dl)  ^(DD) 

q(D1)^(D2)^^^^(DD) 

where 


(ij) 


C.  . » q.  (t  -kT  )q*(t  --t-T  ) 

k4  o a o s 


g(t  -kT  -u)g(t  — t-T  -u)p.  . (u)du 
J o s ’ o s ^13 


(4.43) 


and 


* gd'T.-'tTj  + ^ 


a a N 


-kj^  s k.i  s 


1 a:  i,j  s D 


",  'V"®> 

ni*l 


(4.44) 


The  average  BER  haa  the  aaine  form  aa  4.40 

-1 

G C| 
o 


id  “ 7(det|l  + ^ g“^  c|) 


(4.45) 


but  the  matrix  rank  haa  been  increaaed  from  the  number  of  forward 
filter  weighta  K to  DK  where  D ia  the  number  of  diver aity  channela. 
Fortunately  for  the  uaual  euigle  diveraity  application  the  matrix 
rank  can  be  reduced  to  2K.  Under  the  conditiona  apecifled  in  4.3 
for  a dual  vertical  angle  diveraity  ayatem  in  conjunction  with 
apace  and/or  frequency  diveraity,  the  matricea  C and  G reduce  to 
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G 


(4.46) 


r 


e 


c(ll) 

c(12)  ^(22) 


G 

G 


(11) 

(12) 


.(12)  ' 
g(22) 


and  the  mean  BER  is 

- T 

where  0/2  is  the  ninriber  of  main  beaun  diversity  channels. 

Equation  (4.47)  with  definitions  (4.43)  and  (4.44)  for  the 
matrices;  i^  » 1,  2;  is  the  basic  equation  used  in  this  progreun 
to  compute  mean  BER  for  angle  diversity  performance  of  digital 
systems. 

4.2  Performamce  Prediction  Results 

The  method  of  computing  the  average  bit  error  probability 
for  a troposcatter  equalizer  has  been  applied  to  the  three  tropo- 
scatter  system  examples.  The  path  parameters  for  these  systems 
have  been  given  in  Table  3.1  of  Section  3.  The  troposcatter 
equalizer  chosen  for  the  analysis  is  the  MD-918  modem  developed 
by  Sylvania/SIGNATRON  for  U.S.Army  ECOM  under  the  Megabit  Digital 
Troposcatter  Subsystem  [4.3].  The  modem  uses  differentially  en- 
coded 4 PSK  modulation  and  a four  channel  adaptive  decision- 
feedback  equalizer  [4.4]  for  reception.  There  are  three  taps 
per  forward  filter  diversity  channel  with  tap  spacing  equal  to 
one-half  a 4 PSK  symbol  Interval.  The  backward  filter  uses  three 
taps  in  its  trzuisversal  filter  realization  with  tap  spacing 
equal  to  the  symbol  Interval. 

Before  the  error  probability  can  be  computed,  squint  loss 
must  be  calculated  and  the  29  multipath  spread  must  be  determined 


% -1 

1 + 5-  ® =1; 

o 


(4.47) 
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Bores ight/ 
Horizon  Angle 

Squint 

Loss 

Main 

Beeun 

Elev. 

Becun 

Product 

Be5un 

Beamwidth 

dB 

nsec 

nsec 

nsec 

1/2 

1.9 

132. 

164. 

130. 

1/4 

5.4 

155. 

273. 

200. 

1/4 

4.3 

201. 

329. 

240. 

for  the  main  beam,  elevated  beam,  and  cross  beeun  (or  product) 
multipath  profiles.  These  values  have  been  confuted  in  Section 
3 and  were  presented  in  Figs.  3.2  - 3.4  and  3.8  - 3.10.  The 
selected  values  for  the  error  probability  computation  are  given 
in  Tcdsle  4.1 

Table  4.1 

2a  Multipath  Spread 


System 


jRAOC 

I S.Tepesi/ 

: Yamanlar 

• Oslo/ 

: Kristiansand 


The  data  rate  chosen  for  the  calculation  was  6.3  Mb/s  as  it 

corresponds  to  typical  user  requirements  for  strategic  digital 

troposcatter.  With  a 4 PSK  modulation  format,  the  bit  symbol 

interval  T is  then  317  nanoseconds.  The  2a  /T  values  for  the 

o 

main  beam  diversity  are  0..42,  0.49,  and  0.63  for  the  RADC, 
S.Tepesi,  and  the  Oslo  path,  respectively. 

The  results  of  the  error  probability  confutation  are  plot- 
ted as  a function  of  the  mean  diversity  ^or  no  diversity 

reception  (1  power  anflifier).  For  an  angle  diversity  applica- 
tion, the  squint  loss  and  the  availability  of  a redundant  power 
amplifier  must  be  considered  in  any  system  comparison  with  con- 
ventional systems.  Since  the  2S/2F  configuration  requires  two 
power  amplifiers, any  comparison  with  angle  diversity  should  be 
on  the  same  basis  and  it  also  follows  that  a 2S/2A  system  would 
require  a second  power  suiflifler  for  failure  and  maintenance 
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redundancy  reasons.  For  these  practical  considerations  the  2S/ 

2A  system  utilizes  two  power  amplifiers  tuned  to  the  same  fre- 
quency. Thus,  in  effect,  the  transmit  power  per  diversity  of  a 
2S/2A  system  is  3 dB  more  than  in  analogous  2S/2F  system.  This 
factor  is  not  considered  in  the  presentation  of  error  probability 
vs.  since  represents  received  bit  energy.  The  3 dB 

advantage  of  the  2S/2A  system  over  the  2S/2F  system  is  included 
in  a later  comparison  of  average  bit  error  probability  vs. time 

availability.  For  a 2S/2A  configuration  the  squint  loss,  S_  in 

•p.  I# 

dB,  is  included  in  the  mean  diversity  as  follows; 

Mean  Diversity  * No  Diversity  Ej^/N^  - 

which  accounts  for  all  of  the  elevated  be2uns  having  a loss  of 
S-  dB  relative  to  the  main  beeun. 

Id 

The  use  of  mean  diversity  ^/N^  allows  an  evaluation 
of  the  effects  of  diversity  correlation  and  implicit  diversity 
in  the  compeurison  between  frequency  and  angle  diversity.  The 
correlation  between  diversity  branches  reduces  the  effectiveness 
of  angle  diversity  but  for  small  to  moderate  multipath  condi- 
tions, the  increased  multipath  spread  of  the  elevated  beams 
produces  more  implicit  diversity  in  the  angle  diversity  system. 
Thus,  the  predicted  results  shown  in  Figs.  4.3  - 4.5  show 
slightly  better  performance  for  the  2S/2A  system  when  2a^/T  for 
the  malhbeam  diversity  is  greater  than  0.49  (Figs. 4. 4, 4. 5)  but 
slightly  better  performance  for  the  2S/2F  system  when  the  cor- 
relation effect  is  dominant  at  smaller  multipath  spreads  (Fig. 
4.3).  For  larger  multipath  spread/ syirbol  length  ratios,  it  can 
be  expected  that  the  2S/2F  system  will  again  outperform  the 
2S/2A  system  because  the  multipath  spread  in  the  elevated  beam 
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will  exceed  the  design  capability  of  the  modem  equalizer.  This 
result  is  shown  in  Fig.  4.6  where  predicted  performance  for  a 
12.6  Mb/s  system  operating  on  the  172  mile  L-band  path  is  shown. 
It  should  be  pointed  out,  however,  that  an  expanded,  i.e.,  more 
taps,  DFB  could  be  used  to  compensate  for  this  effect  with  the  re- 
sult of  better  angle  diversity  performance  than  frequency  diver- 
sity performance  for  these  multipath  ranges. 


It  should  be  emphasized  again  that  the  performance  comparisons 
gleaned  from  Figs. 4. 3 - 4.6  are  based  on  an  equal  received  bit 
energy  constraint.  This  comparison  is  useful  in  determining  the 
relative  effects  of  diversity  correlation  and  implicit  diversity. 
It  can  be  very  misleading  if  this  constraint  is  used  for  system 
comparisons.  A realistic  method  of  comparing  system  performance 
is  to  compare  the  average  bit  error  probability  as  a function  of 
time  availability  for  each  system,  i.e.,  what  fraction  of  the 
time  for  a given  system  will  the  average  bit  error  probability 
be  larger  than  a certain  value.  To  accomplish  this  comparison 
we  first  determine  the  distribution  of  the  median  path  loss  using 
the  common  volvune  integration  model  of  Section  2 and  the  NBS 
time  availaibility  analysis  reviewed  in  Section  3.  The  common 
volume  integration  model  includes  the  effect  of  aperture  to 
medium  coupling  loss  and  provides  a path  loss  parameter  normal- 


ized to  unit  gain  antennas.  To  convert  median  path  loss  into 


Ej^/N^  in  dB, 


the  relationship  is 


ij^/N^  » P^+  G^+  L - NF  - lOLOG(Rj^)  + 174  + 1.6 
where 


" transmit  power  in  dBm 

^T/R  * antenna  gain  of  transmitter  receiver  in  dB 
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/ji/e-/?rt(SE  Bit  f^ffcemuTY 


L • total  path  loss  including  coupling  loss  effects 
NF  > noise  figure 

» data  rate  in  b/s 

174  « noise  power  in  1 Hz  bandwidth  in  dBm 

1.6  s dB  amount  that  the  mean  exceeds  the  median  for  a 
Rayleigh  fading  channel. 

The  antenna  gain  formula  for  troposcatter  reflectors  of 
diameter  D is 

G - 10  LOG^TT^D^e/X^) 

where  0 < e < 1 is  the  antenna  efficiency  with  typical  value  on 
the  order  of  0.57. 

In  the  calculation  of  median  path  loss  the  refractive  index 
variance  dependence  is  based  on  Fried's  [4.5]  results  at  optical 
wavelengths.  Also  a service  probability  of  0. 95  has  been  selected 
as  a conservative  measure.  We  have  not  collected  sufficient 
empirical  data  to  establish  the  eibsolute  accuracy  of  the  pre- 
diction model.  However,  the  model  is  very  useful  for  comparisons 
of  different  systems  on  a particular  link.  During  the  next  study 
phase,  effort  is  planned  to  determine  the  absolute  accuracy  of 
the  median  path  loss  prediction  over  a wide  range  of  troposcat- 
ter system  types..  For  the  RAOC  test  link  where  empirical  data 
is  available,  the  winter  median  path  losses  for  the  two  space 
diversity  antenna  receivers  were  measured  [4.6]  to  be  258  and 
260  dB.  Our  predicted  value  for  this  link  for  winter  conditions 
is  258  dB  which  is  excellent  agreement.  The  path  loss  test  re- 
port [4.7]  for  the  link  acceptance  of  the  Oslo-Krlstlansand  link 
in  the  summer  of  1962  shows  a median  path  loss  of  217  dB.  Be- 
cause considerably  stronger  signals  are  experienced  in  the 


4-27 


suimner,  this  result  also  compeures  favorably  with  our  winter 
prediction  of  231  dB  for  this  link.  No  empirical  data  were 
availedile  on  the  S.Tepesi-Yaunanlar  DCS  link. 


The  total  path  loss  (or  distributions  for  the  three 

exaunple  systems  are  given  in  Figs.  4.7  - 4.9.  The  average  bit 
error  probability  curves  from  Figs.  4.3  - 4.5  are  used  with  the 
E|^/N^  distribution  to  determine  the  error  rate  availedsility. 
Note  that  the  ordinates  of  Figs.  4.7  - 4.9  represent  the  no 
diversity  received  bit  energy  with  a single  power  eimplifier 
whereas  the  ^d^scissas  of  Figs.  4.3  - 4.5  are  expressed  as  the 
average  over  the  diversity  branches  of  the  received  bit  energy. 
The  relationships  between  these  values  are  given  by  the  follow- 
ing- for  the  different  system  configurations: 

ij^/N^(lPA)  = AVG  Ej^/N^  2S/2F 

Ij^/N^(1PA)  * S^/2  - 3 + AVG  2S/2A 

Ej^/N^(1PA)  » S^/2  + AVG  Ej^/N^  2S/2F/2A 


The  3 dB  for  2S/2A  accounts  for  the  use  of  the  additional 
power  amplifier  at  the  seume  frequency  in  this  configuration. 

The  quantity  S.  is  the  squint  loss  of  the  elevated  beeun  relative 

id 

to  the  main  beaun.  Since  in  many  examples  the  2S/2F  and  2S/2A 
performance  are  approximately  equal  on  an  AVE  basis,  the 

system  gain  of  2S/2A  over  2S/2F  is  equal  to  3 - S^/2  which  is 
usiially  positive  and  yields  the  system  inqprovements  shown  in 
Figs.  4. 7-4. 9.  Note  that  the  performance  improvement  afforded 
by  2S/2P/2A  over  2S/2A  is  not  all  that  large  because  the  added 
diversity  is  reduced  by  the  need  to  utilize  the  second  power 
an^lifier  at  the  diversity  frequency  resulting  in  a 3 dB  system 
loss. 
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Predicted  System  Performance,  Oslo-Kristiansand 


Most  DCS  troposcatter  systems  have  operating  carrier  fre- 
quencies below  1 GHz.  In  certain  situations,  the  avalleJ^lllty 
of  frequencies  within  the  host  country  might  make  It  attractive 
to  convert  a troposcatter  link  to  an  S-band  or  C-band  carrier 
frequency.  At  these  higher  frequencies  the  path  loss, aperture 
to  medium  coupling  loss,  and  antenna  gains  are  larger.  For  the 

turbulent  scattering  model  the  scatter  angle  dependence  of  the 

-11/3 

scatter  power  varies  as  9 which  results  In  a path  loss  de- 

pendence (not  Including  coupling  loss)  with  frequency  of  . 
This  dependence  should  be  contrasted  with  the  NBS  model  which  Is 
based  on  a layer  theory  with  a path  loss  dependence  of  f^.  The 
turbulent  scatter  theory  Is  viewed  as  more  accurate  for  shorter 

wavelength  systems  and  for  worse  tlme-of-the  year,  l.e.,  winter, 

2 

predictions.  The  antenna  gain  varies  as  f and  when  both  trans- 
mitting are  receiving  antennas  are  taken  Into  account  the  total 

4 

antenna  gain  improvement  Is  f . Thus  In  the  absence  of  coupling 
loss,  higher  frequencies  of  operation  would  result  In  an  Improve- 
ment of  . The  coupling  loss  results  from  a shrinkage  of  the 

usable  common  volume  and  thus  Is  a function  of  the  beamwldth 

-11/3 

and  the  scatter  power  dependence  9 . For  very  narrow  beam- 

width,  n « 9 , where  9 Is  the  smallest  scattering  emgle  visible 
o o 

from  both  terminals,  the  scatter  power  dependence  Is  not  in^rtant 
I and  the  coupling  loss  Is  solely  a function  of  the  beamwldth. 

For  this  asynptotlc  small  beamwldth  case,  the  coupling  loss  Is 
proportional  to  the  comnon  volume  shrinkage  which  varies  as  n" 


f 
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However,  for  L-band  systems,  the  asymptotic  limit  is  not 
approached  and  one  might  realistically  expect  the  link  SNR  to 
improve  with  frequency.  Moreover,  one  must  consider  the  dis- 
tribution of  total  path  loss  (total  path  loss  includes  coupling 
loss  effects)  at  the  different  frequencies  of  interest.  In  many 
examples  the  NBS  model  for  long  term  variability  predicts  smaller 
variation  in  path  loss  for  higher  frequency  systems.  This  smaller 
variation  would  result  in  a better  system  availability  at  higher 
frequency. 

A quantitative  evaluation  of  the  above  results  from  com- 
putation of  total  path  loss  for  the  2 example  L-band  systems 
after  hypothetical  conversion  to  S and  C band.  The  results  are 
given  in  Table  4.2.  For  convenience,  ^ Mb/s  rate 

is  used  to  define  the  link  SNR.  We  see  that  for  both  examples 
conversion  to  S-band  provides  considerable  improvement  which 
is  even  larger  at  the  higher  availabilities.  Conversion 
from  S-band  to  C-band  provides  about  the  same  performance  as 
the  asyn^totic  limit  is  still  not  achieved. 

Although  these  results  are  very  encouraging  for  applica- 
tions where  frequency  conversion  are  being  considered,  there 
are  two  factors  which  the  cautious  engineer  should  consider. 

First,  insufficient  empirical  data  has  been  collected  amd 
analyzed  to  validate  the  turbulent  scatter  model  and  the  result- 
ing f dependence  of  path  loss  over  the  frequency  range  of 
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Table  4.2 

E]|j/Nq  Frequency  Dependence 
6.3  Mb/s,  Service  Prob.  =>  0.95 


Avail- 

ability 

172  Mile,f  » 0.87  GHz 
c 

60'  Ant.,  10  KW 

205  Mile,f  a 0 
c 

66*  Ant.,  10 

.9  GHz 

KW 

^c 

3fc 

5fc 

fc 

3fc 

5fc 

Vc 

V^o 

Vo 

Vo 

90% 

23.5 

27.9 

27.3 

22.5 

24.5 

24.1 

99% 

15.5 

21.5 

20.9 

15.5 

19.1 

18.7 

99.9% 

9.1 

16.6 

16.0 

10.4 

15.0 

14.6 

j99.99% 

t 

4e6 

I2e4 

L__^ 

11.8 

6.2 

11.5 

11.1 

interest.  Second,  the  NBS  long  term  variability  model  is  based 
upon  en^irlcal  data  taken  primarily  on  systems  operating  below 
1 GHz. 
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APPENDIX 

POWER  FADING  CORRELATION  COEFFICIENT 

The  power  fading  correlation  coefficient  can  be  measured 
from  the  received  powers  at  watts  or  in  dBm.  In  this  Appendix 
we  derive  the  relationship  between  these  coefficients.  It  is 
shown  for  a special  case  that  the  watt  measure  of  the  correla- 
tion coefficient  is  overbounded  by  the  dBm  measure  and  for  a 
practical  exan^lethe  coefficients  are  shown  to  be  close  in  value. 

The  density  function  for  the  hourly  median  path  loss  (or 
received  power)  in  dB  (or  dBm)  is  given  by 

1 -x’r“^x/2ct^ 

f(x  X ) = r e (A.l) 

2TTa|Rr 

where 


m^  ■ mean  value  of  hourly  median  path  loss  (dBm)  for  mainbeam 

m^  * mean  value  of  hourly  median  path  loss  (dBm)  for  eleva- 
ted beam 

o ■ standard  deviation  of  mainbeam  path  loss 
aa  " standard  deviation  of  elevated  beam  path  loss 

p ■ power  fading  correlation  coefficient. 


The  correlation  coefficient  is  defined  as 


( *1^  - (*2^  - ^ 


5 -(-1*2  - "■i”2)/»»^ 


(A. 2) 


The  path  loss  in  watts  is 


X./20  0.115  X.  Ox 


y.  * 10 
1 


^ = e ^ , a = 0.115,  i « 1,  2 


and  its  correlation  coefficient  is 


/ 2 - 2 - 2\h  1 

(yi  - Yi  j^(y2  “ ^2  ) 


The  mean  and  variance  for  is  found  by  direct  integration 


(A. 3) 


V.  = 1 1 e 

— *<„ 

oJ-2^ 


2 2 2 2 
ax-(x-m  ) /2a  om^^+isa  a 

dx  " ® 


O *9 

" - , »2,_  2 2am-+2a  a 

2 1 r 2ax-(x-m, ) /2a 1 

y = — — I e 1 - e 

aV2TT  -®  dx 


2 2 2 2 
2ow,  + a a 2a  a 


2 2-2  1 / 

^1  - Yl  - Yl  » e [ 


The  results  for  y^  follow  immediately  by  substituting  aa 
for  a in  the  above  expressions. 


A-2 


T 

r 


_ om^  + 
72  - e 


2 2 2 

- 2am^  + 2a  a a 

2 2 

^2  ■• 


— -2  ^on^+aVo^  2 2^2 

"2  * ^2  - ^2  ■ ® (®  -1) 


To  find  yj^y2  » w«  find  it  convenient  to  perform  an  ortho- 
gonal transformation  on  x to  produce  au  independent  joint  normal 
process.  Let  Q be  the  orthogonal  matrix  vdiich  diagonalizes  R, 

i . e . f 

QW  . r , r.^  . i,  j . 1,  2 

where 


and  e^j^  and  are  the  eigenvectors  of  R, 
Re^  * 1 i””  * 


The  transformation  is 

u - Qx,  X » Q'u  (A. 4) 

and  the  density  function  of  u is 


A-3 


f(u) 


_1 -u'T"^3^2o^ 

\ . € 


(A. 5) 


2tto1r 


which  shows  the  independent  nature  of  and  u^.  Now 


a(x,  + x_  ) 


YlYj  » e 


12'  b'x  + b'm 


where 


b * 


m 


a 

a 

m. 


U “2 


Since  x * Q'u*  we  have 


b'Q'u  b’m  c'u  b'm  *^2^2 

y^2  ®®  ® ® ”®  ® ® 


where  we  have  used  the  independence  of  u^  and  u^  and  where 


Qb 


GtX 


We  have  already  evaluated  terms  of  the  form  e to  find 


•=i”i 

e * e 
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The  power  fading  correlation  coefficient  is  thus 


+a(m  +m2)  J"  ” (<=1  Ha^o^(l+a^)  " 

i Le 2S :i_ 

a(m  +m  ) / «2_2/  2_2  \ „2  2 2/  „2  2 2 N 

> ^ W “ U'"  ^ -1)  a O a O 


-c(e  \ +c  X )/2 
e '' -e 


(A. 6) 


For  the  special  case  of  equal  variances,  a«l,  (A. 6)  reduces 
to  a much  simpler  fom.  The  eigenvalues  and  eigenvector  for  R 
when  a-l  are 


» a‘(i  - p) 


\2  = ^ a + P) 


% “ 


V2 

1 

‘i 


-1 

U.  J 
1 
1 


with  the  result 


2 2 
a o p 


-1 


2 2 
a a 

e -1 


(A. 7) 


A-5 


p < 1 - 

- - 4 4 6 6 

2 2 , a a , 

2 6 

Xl  + (xV  + 1 

p < 1 - 

22  44  <l-e.  , 

1 + J.'  + A_2_  + 

Numerical  examples  show  that  to  a good  approximation 
p * p for  the  above  special  case  and  when  a t !•  For  example, 
for  the  BTL  data  we  have 

p - 0.74 
o • 6.1 
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